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Abstract. Growth estimates of complex orthogonal polynomials with respect to the area measure 
supported by a disjoint union of planar Jordan domains (called, in short, an archipelago) are 
obtained by a combination of methods of potential theory and rational approximation theory. 
The study of the asymptotic behavior of the roots of these polynomials reveals a surprisingly rich 
geometry, which reflects three characteristics: the relative position of an island in the archipelago, 
the analytic continuation picture of the Schwarz function of every individual boundary and the 
singular points of the exterior Green function. By way of explicit example, fine asymptotics are 
obtained for the lemniscate archipelago \z m — 1| < r m , < r < 1, which consists of m islands. The 
asymptotic analysis of the Christoffel functions associated to the same orthogonal polynomials leads 
to a very accurate reconstruction algorithm of the shape of the archipelago, knowing only finitely 
many of its power moments. This work naturally complements a 1969 study by H. Widom of Szego 
orthogonal polynomials on an archipelago and the more recent asymptotic analysis of Bergman 
orthogonal polynomials unveiled by the last two authors and their collaborators. 



Archipelago n. (pi. archipelagos or archipelagoes) an extensive group of islands. 

1. Introduction 

The study of orthogonal polynomials, resurrected recently by many groups of scientists, some 
departing from the classical framework of constructive approximation to fields as far as quantum 
computing or number theory, does not need an introduction. Maybe only our predilection in the 
present work for complex analytic orthogonal polynomials on disconnected open sets needs some 
justification. 

Complex orthogonal polynomials naturally came into focus quite a few decades ago in connection 
with problems in rational approximation theory and conformal mapping. The major result, pro- 
viding strong asymptotics for Bergman orthogonal polynomials in a domain with analytic Jordan 
boundary, goes back to 1923 to a landmark article by T. Carleman [HJ. About the same time S. 
Bernstein discovered that the analogue of Taylor series in non-circular domains (specifically ellipses 
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in his case) is a Fourier expansion in terms of orthogonal polynomials that are well adapted to 
the boundary shape, a phenomenon later elucidated in full generality by J.L. Walsh [36]. Then, it 
came as no surprise that good approximations of conformal mappings of simply-connected planar 
domains bear on the Bergman orthogonal polynomials, that is those with respect to the area mea- 
sure supported by these domains. By contrast, the theory of orthogonal polynomials on the line or 
on the circle has a longer and glorious history, a much wider area of applications and has attracted 
an order of magnitude more attention. For history and details the reader can consult the surveys 
[22] and [33] or the monographs [51 [24l [26| [30] . 

Bergman orthogonal polynomials provide a canonical orthonormal basis in the Bergman space 
of square summable analytic functions associated to a bounded Jordan domain of the complex 
plane. Contrary to the Hardy-Smirnov space, that is roughly speaking the closure of polynomials 
in the L 2 space with respect to the arc-length measure on a rectifiable Jordan curve, the functions 
belonging to the Bergman space do not possess non-tangential values on the boundary. This makes 
their study much more challenging, and less complete as of today. For instance, it is of recent 
date that the analogues of Blaschke products associated to the Hardy space of the disk have been 
discovered: the so-called contractive divisors in the Bergman space of the disk, see the monograph 
by Hedenmalm, Korenblum and Zhu |10j . 

It is our aim to discuss in the present work nth-root and strong estimates for Bergman orthogonal 
polynomials on an archipelago, the asymptotics of their zero distribution, and a reconstruction 
algorithm of the archipelago from a finite set of the associated power moments. The specific choice 
of the above problems and degree of generality were dictated by the present status of the theory of 
complex orthogonal polynomials. 

A brief description of the subjects touched in this article follows. Let G = UjL^Gj be an 
archipelago, that is a finite union of mutually disjoint bounded Jordan domains of the complex 
plane. The Bergman orthonormal polynomials with respect to the area measure supported on G: 

P n (z) = \ n z n + ■■■ , A n >0, 71 = 0,1,2,..., 

carry in a refined (one would be inclined to say, aristocratic) manner the information about G. For 
instance, simple linear algebra provides a constructive bijection between the sequence {P n }%Lo and 
the power moments (correlation matrix entries) 

Mmn(G):= / z n z m dA, m,n>0, (1.1) 
Jg 

where dA stands for the area measure on C. Three major features distinguish Bergman orthogonal 
polynomials: 

(i) An extremality property: P n /X n is the minimum L 2 (G, cL4)-norm monic polynomial of 
degree re, 

(ii) the Bergman kernel K(z, Q = Pj(C)Pj{ z ) collects into a condensed form the (deriva- 
tives of the) conformal mappings from the disk to every connected component Gj, 

(iii) the square root of the Christoffel function A n (z) := 1 / yJ'Y^Jj=Q I Pj ( z ) 1 2 i s the extremum 

value min||g|| L 2 (GiCL4) , q(z) = 1, deg q < re. 

We repeatedly use the above characteristic properties, by combining them with general methods 
of potential theory and function theory. An important object in our work is the multi-valued 
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function 

$0) = exp{#o0,oo) +ign{z)}, zeC\G, 

where gn(z, oo) is the Green function of the exterior domain f2 := C\G, with a pole at infinity, and 
g^ is any harmonic conjugate of gn- We designate the name Walsh function for <3?. At a critical 
moment in our proofs, we rely on the pioneering work of Widom [SB] that refers to Szego's orthogonal 
polynomials on G and their intimate relation to the Walsh function <3?. Our Bergman space setting, 
however, departs in quite a few essential points from the Hardy-Smirnov space scenario. Both 
estimates of the growth of P n (z) and the limiting distribution of the zero sets of {P n }^ =1 depend 
heavily on <3? and its analytic continuation across dG. 

While the estimates for P n (z) are more or less expected, and only how to prove them might bring 
new turns, the zero distribution picture on an archipelago is full of surprises. The uncovering of 
this rich geometry began a few years ago, in the work of two of us and collaborators, on the zero 
distribution of Bergman orthogonal polynomials on specific Jordan domains, cf. [Ill 116} [23] . For 
example, for the single Jordan region consisting of the interior of a regular m-gon, all the zeros of 
P n , n = 1, 2, . . ., lie on the m radial lines joining the center to the vertices, for m = 3 and m = 4 
(sec |13j), while for m > 5 every boundary point of the m-gon attracts zeros of P n , as n — > oo (see 
Thm. 5]). 

As a byproduct of the estimates we have obtained for A n (z), we propose a very accurate 
reconstruction-from-moments algorithm. In general, moment data can be regarded as the arche- 
typal, indirect discrete measurements available to an observer, of a complex structure. To give a 
simple indication how moments appear in geometric tomography, consider a density function p(x, y) 
with compact support in the complex plane. When performing parallel tomography along a fixed 
direction 9, one encounters the values of the Radon transform along the fixed screen 

R(p)(t,9) = (p(x,y),5(xcos9 + ysin9 - t)) 

where 5 stands for Dirac's distribution and (•, •) is the pairing between test functions and distribu- 
tions. Computing then the moments with respect to t yields 

a k (9) = / t k R(p)(t,9)dt = / (xcos9 + y sin 9) k p(x,y)dxdy. 
Jr Jr 2 

Denoting the power moments (with respect to the real variables) by 

°~j,k= x J y k p(x,y)dxdy, i,j>0, 
Jr 2 

we obtain a linear system 

= Yl ( \ ) cos^sin fc -'0 a itk -i. 

After giving 9 a number of distinct values, and noticing that the determinant of the system is 
non-zero, one finds by linear algebra the values Wj,k}^k = Q- This procedure was used by the first 
two authors of this paper in an image reconstruction algorithm based on a different integral trans- 
form of the original measure, see [7] and [8]. In a forthcoming work we plan to compare, both 
computationally and theoretically, these two different reconstruction-from-moments algorithms. 

The paper is organized as follows: Sections [2] and [3] are devoted to necessary background infor- 
mation. We introduce there the notation, conventions and recall certain facts from potential theory 
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and function theory of a complex variable that needed for the rest of the work. Sections 0] (Growth 
Estimates), [5] (Reconstruction of the Archipelago from Moments) and [6] (Asymptotic Behavior of 
Zeros) contain the statements of the main results. In Section[7]we enter into the only computational 
details available among all archipelagoes: disconnected lemniscates with central symmetry. Finally, 
Section [8] contains proofs of previously stated lemmata, propositions and theorems. 

2. Basic concepts 

2.1. General notations and definitions. The unit disk, the exterior disk and the extended 
complex plane are denoted, respectively, 

B := {z E C : \z\ < 1}, A := {z £ C : \z\ > 1} U {oo}, C:=CU{oo}. 

For the area measure in the complex plane we use dA = dA(z) = dxdy, and for the arc-length 
measure on a curve we use \dz\. By a measure in general, we always understand a positive Borel 
measure which is finite on compact sets. The closed support of a measure /i is denoted by supp fj,. 

As to curves in the complex plane, we shall use the following terminology: a Jordan curve is a 
homeomorphic image of the unit circle into C. (Thus, every Jordan curve in the present work will 
be bounded.) An analytic Jordan curve is the image of the unit circle under an analytic function, 
defined and univalent in a neighborhood of the circle. Thus an analytic Jordan curve is by definition 
smooth. We shall sometimes need to discuss also Jordan curves which are only piecewise analytic. 
The appropriate definitions will then be introduced as needed. 

If L is a Jordan curve, we denote by int(L) and ext(L) the bounded and unbounded, respectively, 
components of C \ L. By a Jordan domain we mean the interior of a Jordan curve. If E C C is 
any set, Co (E) denotes its convex hull. 

The set of polynomials of degree at most n is denoted by V n . 

2.2. Bergman spaces and polynomials. The main characters in our story are the Bergman 
orthogonal polynomials associated to an archipelago G := Uj =1 Gj, where G\,...,Gn are Jordan 
domains with mutually disjoint closures. Set Tj := dGj and T := U^LjTj. For later use we introduce 
also the exterior domain $7 := C \ G. Note that T = dG = dQ,. 

Let {P n }^L denote the sequence of Bergman orthogonal polynomials associated with G. This is 
defined as the sequence of polynomials 

P n (z) = \ n z n + ■ ■ ■ , A n >0, n = 0,l,2,..., 

that are obtained by orthonormalizing the sequence 1, z, z 2 , . . . , with respect to the inner product 

</,<?) := / f{z)W)dA. 
Jg 

Equivalently, the corresponding monic polynomials P n (z)/X n , can be defined as the unique monic 
polynomials of minimal L 2 -norm over G: 

\\^Pn\\L\G)=m n {G,dA) := min \\z n + r{z)\\ L 2 [G) , (2.1) 
where ||/|| i2(G) := (Z,/) 1 / 2 . Thus, 

— = m n (G,dA). 



BERGMAN POLYNOMIALS ON AN ARCHIPELAGO 5 

Let Lq(G) denote the Bergman space associated with G and (•,•): 

L l(G) ■= {/ analytic on G and ||/|| L 2 (G) < 00} . 

Note that L 2 a {G) is a Hilbert space that possesses a reproducing kernel which we denote by K(z, Q. 
That is, K(z, Q is the unique function K(z, Q :GxG^C such that, for all ( G G, K(-, () G L 2 a (G) 
and 

/(C) = (/,^(.,C)>, V/GL2(G). (2.2) 
Furthermore, due to the reproducing property and the completeness of polynomials in L^(G) (see 
Lemma 13.31 below) . the kernel K(z,£) is given in terms of the Bergman polynomials by 

00 

j=0 

We single out the square root of the inverse of the diagonal of the reproducing kernel of G 

y/K{z,z) 

and the finite sections of K(z, Q and A(z): 

n 1 

K n (z,Q :=J2 P i(0Pj(z), K(z) := ■ (2.3) 

j=0 \/ K n{Z,Z) 

We note that the A n (z)'s are square roots of the so-called Christoffel functions of G. 

2.3. Potential theoretic preliminaries. Let Q be a polynomial of degree n with zeros z\ , Z2, ■ ■ ■ ,z n . 
The normalized counting measure of the zeros of Q is defined by 



1 n 

: =-E^' ( 2 - 4 ) 



n 

k=l 

where 5 Z denotes the unit point mass at the point z. In other words, for any subset A of C, 

. . . number of zeros of Q in A 

vq( a ) = • 

n 

Next, given a sequence {cr n } of Borel measures, we say that {<r n } converges in the weak* sense 
to a measure a, symbolically a n — — > cr, if 

/dcr n > y /dcr, n -> 00, 

for every function / continuous on C. 

For any finite positive Borel measure a of compact support in C, we define its logarithmic 'potential 

by 

U%z) := J log^l-da(i). 

In particular, if Q n is a monic polynomial of degree n, then 

= --log|Q n (z)|. 
n 
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Let E C C be a compact set. Then there is a smallest number 7 G R U {+00} such that there 
exists a probability measure /is on E with < 7 in C. The (logarithmic) capacity of £ is defined 
as cap (E) = e~ 7 (interpreted as zero if 7 = +00). If cap (E) > 0, then /i£ is unique and is called 
the equilibrium measure o/E. For the definition of capacity of more general sets than compact sets 
see, e.g., [20] and [24] • A property that holds everywhere, except on a set of capacity zero, is said 
to hold quasi everywhere (q.e.). For example, it is known that U^ 12 = 7, q.e. on E. 

Let W denote the unbounded component of C \ E. It is known that supp(^s) C dW, /i£ = \iqw 
and cap (E) = cap(cW). If cap (E) > 0, then the equilibrium potential is related to the Green 
function gw{ z -> °o) of W, with pole at infinity, by 

U^(z) =log^— -g w (z, 00), z€W. (2.5) 
cap^Zj j 

In our applications dW will be a finite disjoint union of mutually exterior Jordan curves (typically 
E = G or E = r, W = f2, dW = T = <9E, in the notations of Subsection 12. 2p . Then, every point of 
dW is regular for the Dirichlet problem in W [201 Thm 4.2.2] and therefore: 

(i) 

supp/is = dW, (2.6) 

(ii) 

U^(z) = log— , zeS. (2.7) 
cap(E) 

If ji is a measure on a compact set E with cap (E) > 0, the balayage (or "swept measure" ) of /i 
onto <9E is the unique measure v on <9E having the same exterior potential as /i, i.e., satisfying 

U v = W 1 in C \ E. (2.8) 

The potential U v of v can be constructed as the smallest superharmonic function in C satisfying 
U u > in C \ E. Since itself is competing it follows that, in addition to (plgft . U v < in all 
C. 

2.4. The Green function and its level curves. Returning to the archipelago, let ga(z, 00) 
denote the Green function of Q, = C \ G with pole at infinity. That is, <7n(z,oo) is harmonic in 
\ {00}, vanishes on the boundary r of G and near 00 satisfies 

1 1 

50(2,00) = log £ +log 7F - + 0(- rr ), 2 ->oo. (2.9) 

cap(I ) \z\ 

We consider next what we call the Walsh function associated with £1. This is defined as the 
exponential of the complex Green function, 

$(z) := exp{g n (z,oo) + ig^(z, 00)}, (2.10) 

where g^iz, 00) is a (locally) harmonic conjugate of gn(z, 00) in $7. In the single-component case 
N = 1, (|2.10p defines a conformal mapping from O onto A. In the multiple-component case N > 2, 
$ is a multi-valued analytic function in However, |$(z;)| is single-valued. We refer to Walsh [36} 
§4.1] and Widom |38|. § 4] for comprehensive accounts of the properties of We note in particular 
that $ is single-valued near infinity and, since g^ is unique apart from a constant, that it can be 
chosen so that $ has near infinity a Laurent series expansion of the form 

T / n 1 Ql /x 
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We also note that $'(z)/$>(z) = 2dgn(-, oo)/dz is single-valued and analytic in Q, with periods 

1 f 1 f dg n (z,oo) 

Here Tj is oriented as the boundary of Gj and the normal derivative is directed into f2. If r,- is 
not smooth the path of integration in (|2.12p is understood to be moved slightly into 0. Note that 

Next we consider for R > 1 the level curves (or equipotential loci) of the Green function, 

L R ■= {z £ n : gn(z,oo) = log R} = {z e : \$(z)\ = R} (2.13) 

and the open sets 

tt R := {z £ ft : gn(z,cc) > log R} = {z £ Q : \$(z)\ > R} = ext(L R ), 
Q R := C \ tt R = int(L R ). 

Note that L\ = T, i\ = Vt, Q\ = G. It follows from the maximum principle that Q R is always 
connected. The Green function for £l R is given by 

9n R {z,oo) = g n (z,oo) - log R, (2.14) 

hence the capacity of L R (or Q R ) is 

cap(Lii) = i?cap(r). (2.15) 

Unless stated to the contrary, we hereafter assume that N >2, i.e. G consists of more than one 
island. For small values of R > 1, Q R consists of N components, each of which contains exactly 
one component of G, while for large values of R, Q R is connected (with G C Qr). Consequently, 
we introduce the following sets and numbers: 

Qj,R '■= the component of Q R that contains Gj, j = 1, 2, . . . , N. 
Lj lR :=dg j>R , j = 1,2,..., N. 

Rj := sup{i? : Gj,R contains no other island than Gj}. 

R' := min{i?i, . . . , Rn} = sup{i? : Qr has N exactly components }. 

R" := infjii : Qr is connected } = m.{{R : Q,r is simply connected }. 

Thus, when 1 < R < R' , Qr is the disjoint union of the domains Qj t R, j = 1, 2, ... ,N and Lr 
consists of the ./V mutually exterior analytic Jordan curves Lj_R, j = 1,2, . . . , N , while for R > R" , 
we have Q\ tR = Q2,r = • • • = Qn,r = Qr and L R is a single analytic curve. 

It is well-known that gu(z, oo) has exactly N—l critical points (multiplicities counted), i.e., points 
where the gradient V^q(z, oo), or equivalently = 2dgn(-, oo)/dz, vanishes. These critical 

points show up as singularities of L R , which are points of self-intersection. Such singularities must 
appear when L R changes topology. It follows that there are no critical points in Q R i \ G, at least 
one critical point on each L Rj , j = 1,2, ... ,N (one of them is L R >), at least one on L R u and no 
critical point in Q R ». Any Lj R that does not contain a critical point is an analytic Jordan curve. 
In particular, this applies whenever 1 < R < R' or R" < R < oo. 

When R > R" , $ is the unique conformal map of Qr onto Ar := {w : \w\ > R} that satisfies 
(12. lip near infinity. 

In Figure [T] we illustrate the three different types of level curves Lri, Lrh and Lr with R > R" , 
introduced above. 
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Figure 1. Green level curves 



Remark 2.1. The level curves in Figure [T] were computed by means of Trefethen's MATLAB code 
manydisks.m |34j . This code provides an approximation to the Green function g^(z, oo) in cases 
when G consists of a finite number of disks, realizing an algorithm given in [5]. 

Consider now the N Hilbert spaces L^(Gj) defined by the components Gj, j = 1, 2, . . . , N, and 
let K ' (z, C), j = 1, 2, . . . , N, denote their respective reproducing kernels. Then, it is easy to verify 
that the kernel K (z, () is related to K G i (z, () as follows: 



K G i(z,() if z, (€G V 

if z £ Gj, C € G k , j / k. 



(2.16) 



In view of (|2.16p . we can express K(z,Q in terms of conformal mappings <pj : Gj — ► D, j = 
1,2, . . . , N . This will help us to determine the singularities of K(-, Q) and, in particular, whether 
or not this kernel has a singularity on dGj. This is so because, as it is well-known (see e.g. [H p. 
33]), 



K G J(z,() 



¥>'(*M(C) 



7T 



1 - ipj(z)<Pj(C) 



z,(eGj,j = 1,2, 



(2.17) 



By saying that a function analytic in Gj has a singularity on dGj, we mean that there is no open 
neighborhood of Gj in which the function has an analytic continuation. 



3. Preliminaries 

3.1. The Schwarz function of an analytic curve and extension of harmonic functions. 

Let r be a Jordan curve. Then T is analytic if and only if there exists an analytic function S(z), 
the Schwarz function of T, defined in a full neighborhood of T and satisfying 



S(z) = z for z G T; 
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see [3] and [25]. The map z \— > S(z) is then the anticonformal reflection in T, which is an involution 
(i.e., is its own inverse) on a suitably defined neighborhood of T. In particular, S'(z) ^ in such a 
neighborhood. 

If u is a harmonic function defined at one side of an analytic Jordan curve T and u has boundary 
values zero on T, then u extends as a harmonic function across T by reflection. In terms of the 
Schwarz function S(z) of V the extension is given by 

u(z) = -u(S(zj) (3.1) 
for z on the other side of T (and close to T). Conversely we have the following: 

Lemma 3.1. Let T be a Jordan curve and let u be a (real-valued) harmonic function defined in a 
domain D containing T such that, for some constant c > 0, there holds: 

(i) u = on T, 

(ii) \u\ — > c as z — > dD, 
(hi) Vn / in D, 

where V« denotes the gradient of u. Then T is an analytic curve, the Schwarz function S(z) of T 
is defined in all D, and z i— > S(z) maps D onto itself. Moreover, u and S(z) are related by 113.1]) . 
In particular z i— > S(z) maps a level line u = a of u onto the level line u = —a. 

We note that the Schwarz function is uniquely determined by T, but u is not; there are many 
different harmonic functions that vanish on T. A domain which is mapped into itself by z i— ► S(z) 
will be called a domain of involution for the Schwarz reflection. 

Example 3.1. Assume that, under our main assumptions, one of the components of T, say T\, is 
analytic. Then the Green function u(z) = gn(z, oo) extends harmonically, by the Schwarz reflection 
(|3.ip . from O into G\. We keep the notation gn(z, oo) for so extended Green function. Recall that 
the level lines reflect to level lines, so that for R > 1 close enough to one, L^r is reflected to the 
level line 

L x i ={z G Gi : gn{z,oo) = -\ogR} = {zeG 1 : \${z)\ = -} 

' R Jri 
of the extended Green function (and extended Generally speaking, whenever applicable we 
shall keep notations like Lj tP , Gj,p, L p , Q, p etc. for values p < 1 in case of analytic boundaries. 

As was previously remarked, u(z) = g^(z, oo) has no critical points in the region Qi^^ \G±. It 
follows then from (|3.ip that if the Green function extends harmonically to a region G\ \ Q\ )P with 
■g- < p < 1, then it has no critical points there, and the region D = Q x i \ Q\ )P is symmetric with 
respect to Schwarz reflection and is a region of the kind D discussed in Lemma 13.11 

3.2. Regular measures. The class Reg of measures of orthogonality was introduced by Stahl 
and Totik \17\ Definition 3.1.2] and shown to have many desirable properties. Roughly speaking, 
p G Reg means that in an "n-th root sense" , the sup-norm on the support of p and the L 2 -norm 
generated by p have the same asymptotic behavior (as n — ► oo) for any sequence of polynomials of 
respective degrees n. It is easy to see that area measure enjoys this property. 

Lemma 3.2. The area measure dA\c on G belongs to the class Reg. 

Lemma 13.21 yields the following n-th root asymptotic behavior for the Bergman polynomials P n 
in £1. 
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Proposition 3.1. The following assertions hold: 
(a) 

lim \]l n = (3.2) 

n->oo cap(l ) 

(b) For every z S C \ Co(G) and for any z G Co(G) \ G not a limit 'point of zeros of the P n 's, 
we have 

lim \P n (z)\ 1/n = |$0)|. (3.3) 

n— >oo 

The convergence is uniform on compact subsets ofC\ Co(G). 

(c) 

limsup|P n (z)| 1/n = |$(z)|, zeli, (3.4) 



locally uniformly in Q. 

(d) 



IP'(z) &(z) 
lim -^fl = ^fl, (3.5) 



locally uniformly in C \ Co(G) . 

The first three parts of the proposition follow from Theorems 3.1.1, 3.2.3 of [27] and from 
Theorem III. 4. 7 of |24j . in combination with the results of pQ, because is regular with respect to 
the Dirichlet problem; see e.g. [20l p. 92]. The last assertion (d) is immediate from (b). 

Another fundamental property of Bergman polynomials, whose proof involves a simple extension 
of the simply-connected case treated in Theorem 1, Section 1.3 of Gaier [6] is the following. 

Lemma 3.3. Polynomials are dense in the Hilbert space L^(G). Consequently, for fixed ( G G, 

oo 

K(z,Q = ^Pj£)P n {z), (3.6) 

n=0 

locally uniformly with respect to z in G. 

The analytic continuation properties of K(z, play an essential role in the analysis. The fol- 
lowing notation will be useful in this regard. If / is an analytic function in G, we define 

p(f) := sup{-R : / has an analytic continuation to Qr\ . (3-7) 

Note that 1 < p{f) < oo. The following important lemma, which is an analogue of the Cauchy- 
Hadamard formula, is due to Walsh. 



Lemma 3.4. Let f E L A a {G) . Then, 



hmsupK/,^)! 1 /^— -. (3.8) 

n^oo PU) 



Moreover, 

oo 

f(z) = ^2(f,P n )P n (z), 
n=0 

locally uniformly in G p (f) ■ 
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The result is given in Walsh [36} pp. 130-131] (see also |18} Thm 2.1]) for a single Jordan region 
and, as Walsh asserts, is immediately extendable to several Jordan regions. 

By applying Lemma l3.4l to / = K(-, £), and by using the reproducing property (|2.2p . in conjunc- 
tion with (|2.16p and (13.6|) . we obtain: 

Corollary 3.1. Let j be fixed, 1 < j < N. Then for any C £ Gj, 

limsup|P n (C)| 1/n = * Tys = . r / tsG- t 7V\ p T ' ( 3 - 9 ) 

where (as previously defined) Rj > 1 is the largest R such that the component Qj,R of Qr containing 
Gj contains no other G^. In particular, 

limsup|P n (C)| 1/n = l (3.10) 

n— *oo 

if and only if K Gj (•, Co) has a singularity on dGj, for some (and then for every) Co £ Gj. 

The last statement is based on the observation, from (|2.17j) . that the property of K J Co) having 
a singularity on dGj is independent of the choice of Co (within Gj). We remark also that the 
appearance of Rj in (13. 9p is essential since, for R > Rj, the component Qj,R contains an open set 
where K(-, C) is identically zero (recall (|2.16j) ) and hence not an analytic continuation of K Gj (-, C). 
Corollary 13.11 will be further elaborated in Theorem I6.ll 

Corollary 13.11 describes a basic relationship between the orthogonal polynomials {P n (C)}^ =0 and 
the kernel function K(-, C) which will play an essential role in deriving our zero distribution results 
for the sequence {P n }^ =1 . 

Remark 3.1. A well-known result by Fejer asserts that the zeros of orthogonal polynomials with 
respect to a compactly supported measure a are contained in the closed convex hull of supp a. This 
result was refined by Saff [22] to the interior of the convex hull of supper, provided this convex hull 
is not a line segment. Consequently, we see that all the zeros of the sequence {Pn}^^ are contained 
in the interior of convex hull of G. This fact should be coupled with a result of Widom [37] to the 
effect that, on any compact subset E of f2 and for any n G N, the number of zeros of P n on E is 
bounded independently of n. 

3.3. Carleman estimates. We continue this section by recalling certain results due to T. Carle- 
man and P.K. Suetin, regarding the asymptotic behavior of the Bergman polynomials in the case 
where G consists of a single component (i.e. for N = 1). In this case the Walsh function (|2.10p 
coincides with the unique conformal map $ : — > A satisfying (|2. 1 1 j) . 

The first result requires the boundary T to be analytic (hence the conformal map $ has an 
analytic and univalent continuation across T inside G) and is due to Carleman [3]; see also [6[ 
p. 12]. 

Theorem 3.1. Assume that T is an analytic Jordan curve and let p, < p < 1, be the smallest 
index for which $ is conformal in £l p . Then, 



vr cap(r) n+1 

and 



^ = \/ n !* * J l + 0{ P 2n )}, (3.11) 



\l ^&(z)<S> n (z){l + A n (z)}, (3.12) 
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where 

a (a _ J 0(^/E) P n , if *e ft, ( v 

An ^ j " \ 0(l/v^) (p/r) n , if 2 € L n p < r < 1. l ^ ^ j 

The second result which is due to Suetin [301 Thms 1.1 & 1.2], requires that F can be parameter- 
ized with respect to the arc-length, so that the defining function has a p-th. order derivative (where 
p is a positive integer) in a Holder class of order a. We express this by saying F is C p+a -smooth. 
(In particular, this implies that F can have no corners.) 

Theorem 3.2. Assume that F is C( p+1 ) +Q -smooth, with p + a > 1/2. Then, 



and 



\ n = J 1 l±l _L^{i + o(_L_)} (3.14) 

V vr cap(r)"+ l1 y n 2 P +2aJJ> V J 

^) = y^*'(^ n W{l + o(^)}, zeU. (3.15) 

We emphasize that the above two theorems concern only the case when N = 1. We also 
remark that for the case when T is analytic, E. Miha-Diaz [15] has recently derived an improved 
version of Carleman's theorem for the special case when L p is a piecewise analytic Jordan curve 
without cusps. 

3.4. Comparison of area and line integrals of polynomials. The following observation is due 
to Suetin [29]; see also [301 P- 38]. 

Lemma 3.5. Let G be a Jordan domain with C 1+a -smooth boundary. Then there exists a positive 
constant C , with the property that, for every polynomial Q n of degree at most n, there holds 



||Qn||L 2 (r) < C Vn+1 ||Qn||z,2 (G ), 

where \\ ■ ||^2(p) denotes the I? -norm on T with respect to \dz\. 

The proof in [29] uses the following analogue of Bernstein's inequality (a result Suetin attributes 
to S. Yu. Al'per): 

HQnlU 2 (G) <Cn WQnWh^iG) 
and leads to similar inequalities for L p , 1 < p < oo, or uniform norms. 

4. Growth Estimates 

The main results of this article are stated in this and the next three sections. Their proofs are 
given in Section [51 

We recall the notation and definitions in Section [21 in particular the definition of the archipelago 
G := UjL^j consisting of the union of N Jordan domains in C, with boundaries Tj := dGj. We 
also recall that O := C \ G and note T := uf =1 Tj = dG = d£l. 

Theorem 4.1. Assume that every curve Tj constituting F is C 2+a -smooth. Then there exists a 
positive constant C\ such that 



in + 1 
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In addition, if every Tj is analytic, j = 1,2, ... ,N, then there exists a positive constant C2 such 
that 



The following estimate for the diagonal K(z,z), z G G, of the reproducing kernel follows from 
classical estimates for the boundary growth of the Bergman kernel of a simply-connected domain, 
obtained via conformal mapping techniques. More precisely, by using the results for the hyperbolic 
metric presented by Hayman in [9l pp. 682-692], which require no smoothness for the boundary 
curve, and recalling (I2.17|) . it is easy to verify the following double inequality, holding for every j, 
j = l,2,...,N: 

— ^- < K G Hz,z)< =7 -, zeGj, (4.3) 

levrdist 2 ^,^) ~ K ' vrdist 2 ^,^) 3 V ; 

Thus K{z,z) = K G i(z,z), z G Gj, inherits the same estimates and, clearly, the function A(z) := 
1 / a/ K(z, z) satisfies 

0rdist(z,r j ) < A(z) < 40Fdist(z, Tj), z G Gj. (4.4) 

(Above and throughout this article dist(z,Tj) stands for the Euclidean distance of z from the 
boundary Tj.) 

It is always useful to recall that the monic orthogonal polynomials P n (z)/X n , n = 0,1, . . ., satisfy 
a minimum distance condition with respect to the L 2 -norm on G, in the sense that 

I p 

V" = IIt^IU^G) = min \\z n + a n -iz n ~ 1 + ... + a \\ L 2 {G) . (4.5) 

A n A n v ' 00,01,... ,o n _i v ' 

Similarly, the square root of the Christoffel functions A n (z), n = 0, 1, . . ., defined by by (|2.3p . can 
be described as 

A n (z) = min ||p|| L 2 (G) ; (4.6) 

p(«)=i 

cf. Lemma 18.11 below. Based on the above two simple extremal properties, we derive the follow- 
er / ~G~ 

ing comparison between A n (z) and the functions A n 3 (z) := 1/y K n 3 (z, z) associated with each 
individual island Gj. 

Theorem 4.2. For every j = 1,2, . . . , N and any n G N, 

An 3 \z)< A n (z), z€C. (4.7) 
In addition, ifTj is analytic, then there exist a sequence {7n}^Ln> w ^h < j n < 1 and lim 7 n = 

n — >oo 

geometrically, and a number m G N, m > 1, such that for any n G N, 

l -^A mn {z)< A%(z), zeGj. (4.8) 

Let <3?j denote the normalized, like (12. lip , exterior conformal map $j : C \ Gj — > A. The growth 

G 

of A n 3 (z) inside the island Gj is described as follows. 
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Theorem 4.3. Fix j, j = 1,2, ... ,N, and assume that Tj is analytic. Then there exist positive 
constants C\ , C2 and p < 1 such that for any n S N, 

0< A^(*)-A G '(z)< Cil^^rfdist^r^ + ^V z E~Gj\ Qj iP . (4.9) 

Moreover, 



lim nk%(z) = -^- (4.10) 

n^oo \®j(Z) I 

uniformly for z £Tj. 

Furthermore, if every curve constituting T is analytic then 

C 3 < nk n {z) <C 4 , z £ T, (4.11) 

and 

C 5 dist(z, r) 8{z) < k n (z) < n G N, z^G, (4.12) 

yjn\®(z)\ n 

where C?,,C/±,C<o,C§ are positive constants and 

mz)\ 2 -i 1 



S(z) 



y/(n + l)\<S>(z)\^(\<S>(z)\i-l) + l 



An estimate for k n (z) on T which is finer than (|4.1ip . in the sense that it coincides with (|4.10p 
for the case of a single island, and under weaker smoothness conditions on T, is presented in [32], 
where asymptotics of Christoffel functions defined by more general measures on C are considered. 

Finally, we derive the following exterior estimates for Bergman polynomials. 

Theorem 4.4. Assume that every curve constituting T is analytic. Then the following hold: 

(i) There exists a positive constant C, so that 

\Pn(z)\< d]st C {zX) ^\HzT, z^G. (4.13) 

(ii) For every e > there exist a constant C e > 0, such that 

\P n {z)\ > C e Vn|*(.z)r, dist(z, Co(G)) > e. 

Note that in the region Co(G) \ G the orthogonal polynomials may have zeros (as the case of the 
lemniscates considered in Section [7] illustrates) . 



5. Reconstruction of the archipelago from moments 

The present section contains a brief description of a shape reconstruction algorithm. This al- 
gorithm is motivated by the estimates established in the previous sections. The comparison of 
the speed of convergence and accuracy of this approximation scheme with other known ones (see 
e.g. [8]) will be analyzed in a separate work. 

The algorithm is based on the following observations: 



Remark 5.1. 
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(i) From (j4.4j) we see that the function A(z) is bounded from below and above in G by constants 
times the distance of z to the boundary. Consequently, its truncation 

AJz) = 1 = (5.1) 

VELol^WI 2 

approximates the distance function to T in G. Furthermore, on T and in £1 A n decays to 
zero at certain rates, as n — ► oo. More precisely, a close inspection of the inequalities in 
Theorems 14.21 and 14.31 in conjunction with (|4.4p . reveals the following asymptotic behavior 
of A n (z) in C: 

(a) v^dist^T) < A n {z), z £ G; 

(b) A n (z) < Cdist(z, T), z e Gf]Qp, for some < p < 1 and C > ^7r; 

1 

(c) A n (z) x -, zeT; 

n 

(d) A n (z) x , zen. 

yjn\®{z)\ n 

(ii) In order to construct A n we need to have available the finite section {Pq, P±, . . . , P n } of 
Bergman polynomials, and this can be determined by means of the Gram-Schmidt process, 
requiring only the power moments (jl.ip . of degree less or equal than n in each variable. 

(iii) For any n = 1, 2, . . ., all the zeros of P n {z) lie in the interior of the convex hull of G; see 
Remark 13.11 

The expression ixB means that C\B < A < C2B for positive constants C± and C<i. 

Consequently, Remark 15.11 supports the following algorithm for reconstructing the archipelago 
G, by using a given finite set of the associated power moments 



IHj := (z % ,z?) = / zVcL4(z), i, j = 0, 1, . . . ,n. 
JG 

Reconstruction Algorithm 

1. Use an Arnoldi version of the Gram-Schmidt process, in the way indicated in [28J, to con- 
struct the Bergman polynomials {-Pfc}^ =0 from fiij, i,j = 0, 1, . . . ,n. This involves at the k-step 
the orthonormalization of the set {Pq, P\, . . . , Pk-i, zPk-i}, rather than the set of monomials 
{1, z, . . . , z k }, as in the standard Gram-Schmidt process. 

2. Plot the zeros of P n , 11 = 1,2, ... ,n. 

3. Form A n (z) as in \5.1\) . 

4. Plot the level curves of the function A n (x + iy) on a suitable rectangular frame for (x,y) that 
surrounds the plotted zero set. 

Regarding the stability of the Gram-Schmidt process in the Reconstruction Algorithm, we note 
a fact that pointed out in [28]. That is the Arnoldi version of the Gram-Schmidt does not suffer 
from the severe ill-conditioning associated with its ordinary use; see, for instance, the theoretical 
and numerical evidence reported in [19] . This feature of the Arnoldi Gram-Schmidt has enabled us 
to compute accurately Bergman polynomials for degrees as high as 160. We also note that the use 
of orthogonal polynomials in a reconstruction-from-moments algorithm, was first employed in [28J. 
However, the algorithm of [28] is only suitable for the single island case N = 1. 
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Applications of the Reconstruction Algorithm are illustrated in the following six examples. In 
each example, the only information used from the associated archipelago G was its power moments. 
The resulting plots indicate a remarkable fitting, even in the case of non-smooth boundaries, for 
which our theory, as stated in Section HJ does not apply. 




Figure 2. Level curves of Ai o(z + iy), on {(x, y) : -2 < x < 5, -2 < y < 2}, with 
G as in Example 16.11 




Figure 3. Level curves of Ai o(x + iy), on {(x, y) : -4 < x < 4, -2 < y < 2}, with 
G as in Example 16.21 




Figure 4. Level curves of Aioo(x + iy), on {(x,y) 
G as in Example 16.31 case (i) ■ 



-2 < x < 8, 



-2 < y < 2}, with 
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FIGURE 5. Level lines of A 100 (:e + iy), on {(x, y) : -1 < x < 4, -2 < y < 2}, with 
G as in Example 16.51 




Figure 6. Level lines of Ai o(£ + iy), on {(x,y) : -1 < x < 6, — 2 < y < 2}, with 
G as in Example I6.4L 
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Figure 7. Level lines of A n (x+iy), for the values of n (from left two right) 25, 50, 75 
and 100, on {(x,y) : — 1 < x < 4, — 2 < y < 3}, with G formed by the three disjoint 
disks of Example I6.6L 
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6. Asymptotic behavior of zeros 

6.1. General statements. The first result of this section is our general theorem on the asymptotic 
behavior of the zeros of the Bergman polynomials {Pn}^^ on an archipelago of N Jordan domains. 
It is established under the general assumptions made at the beginning of Section 12.21 In particular 
we note that, unlike the theory presented in Section HJ no extra smoothness is required for the 
boundary curves Tj here. The result below, which is valid for any N > 1, requires some special 
attention for the single island case N = 1. 

Theorem 6.1. Consider the following extension of the Green function gn(^oo) to all C: 



(recall j3.7\)) and define 



U(.,oc) zen, 

\-\0g P (K(;Z)), ZGG, 

P = Pg:= t^A/j, (6.2) 

where the Laplacian is taken in the sense of distributions. Let C denote the set of weak* cluster 
points of the counting measures {vp n }^ =1 , i.e., the set of measures a for which there exists a 
subsequence M a C N such that vp n — — > a, as n — > oo, n G J\f a - The following assertions hold. 

(i) The function h is harmonic in ; subharmonic in all C; hence (3 is a positive unit measure 
with support contained in G. In addition, if N > 2, then h is continuous and bounded from 
below. If N = 1, then h can take the value —oo at most at two points, and outside these 
points h is continuous. 

(ii) 

U?{z) = log — ^— - h(z), z G C. (6.3) 
cap (1 ) 

and balayage of (3 onto V gives the equilibrium measure fir of T : 

tn»>u«inc 

[UP = infl. v ' 

("i) i 

limsup — log |P n (^)| = h(z), z G C, (6-5) 

n^oo Tl 

liminf U VB ^(z) = U p (z), z G C. (6.6) 

n— >oo 

Moreover, in C \ Co(G) these equalities hold with limsup and liminf replaced by lim. 

(iv) The set of cluster points C is nonempty, and for any a G C, 

I"'-"" (6.7) 
I U a = U p in the unbounded component of C \ supp (3. 

(v) The measure (3 is the lower envelope of C in the sense that 

C/ /3 = lsc(inf U a ), 

where "Isc" denotes lower semicontinuous regularization. (This means that U@ is the supre- 
mum of all lower semicontinuous functions that are < info-gc U a .) In addition, ifV is any 
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component of C \ supp/3, then for any a G C either U c > U 13 in V or U a = U^ 3 in V; and 
there exists a a G C such that the latter holds. 

(vi) If C has only one element, then this is (3 and 

u Pn — U p, n -> oo, n G N, (6.8) 

i.e., the full sequence converges to [3. 

(vii) Assume that (3 satisfies 

(a) supp/3 is a nullset with respect to area measure, 

(b) C \ supp (3 is connected. 

Then (3 is the unique element inC; hence A6.8\) holds. If (a) holds and (in place of (b)) 

(c) C \ supp (3 has at most two components, 
then f3 G C. 

Remark 6.1. The measure [3 = (3q is canonically associated to G via the Bergman kernel. Con- 
structive formulas for [3g (or rather its potential) will be given in the proof (e.g. (|8.20|) - (|8.22|) ) and 
will be further elaborated in the examples of Section [ 



Remark 6.2. Well-known properties for any a G C follow immediately from (ii) and (iv): That is, 
U a = U^ F in f2, supper C G and balayage of a onto T gives the equilibrium distribution (see 
e.g. jM Thm III.4.7]). 

Remark 6.3. We know of no example where (3 isn't itself in C. However it remains an open question 
whether it is always so. 

Remark 6.4. A measure (3 satisfying (16. 4p together with (a) and (b) in (vii) may be viewed as a 
potential theoretic skeleton for /ir (or "Madonna body" , in view of a common shape of supp (3; cf. 

mills])- 

Remark 6.5. When N = 1, h may take the value — oo at one or two points. Note that, by (6.1), 
h{a) = —oo if and only if K{z,a) is an entire function of z. With G = ID we have h(z) = log \z\, 
i.e., one pole for h. An example with two poles is the following. 

Choose a number A > 1 and let G be the image of the unit disk under the conformal map 

the branch chosen so that "0(0) = 0. The inverse map is 

ip{z) = Atanhz, 

which is meromorphic in the entire complex plane. Here ip maps the disk \w\ < A onto the strip 
|Imz| < j. Hence G, which is the image of \w\ < 1, is a subdomain of that strip (a kind of an 
oval) . 

The function (p does not attain the values ±A anywhere in the complex plane and the set 
ipl^ 1 (l/<f(C)), which will play an important role in the proof of the theorem, may therefore be 
empty for up to two values of C G G. In fact, this occurs for £ = ±a G G, where a = \ log jcs^i > 0- 

At these points, h(±a) = — oo, K(z, ±a) = - ~ 1 e ±2z . One also finds that (3 is a measure supported 
on the line segment [—a, a] and hence is a Madonna body. 
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We call a boundary curve Tj singular if some conformal map ipj : Gj — > B does not extend 
analytically to a full neighborhood of Gj, i.e., if p(tpj) = 1, or equivalently if p(K(-, z)) = 1, z G G^; 
see ()2.16p and (|2.17p . Clearly, this property is independent of the choice of the conformal map 
(fj. Note that a boundary component that is not singular in the above sense still need not be 
fully smooth: it may be piecewise analytic but have certain kinds of corners so that (fj extends 
analytically across Tj but the extension is not univalent. This would be the case, for instance, if 
Gj is a rectangle. 

Corollary 6.1. For each j = 1, . . . , N, the following statements are equivalent: 

(i) Tj is singular. 

(ii) /%. =Pt\g 3 - 

(iii) There is a subsequence M = Afj C N such that, with V any neighborhood of Gj not meeting 
the other islands (e.g., V = Qj^), 

vp n \v Pt\vi n — ► oo, n£j\f. (6.9) 

Clearly, under the conditions of the above corollary a certain proportion of the zeros of the 
Bergman polynomials converge to the part of the equilibrium measure located on Tj. By a reasoning 
as in deriving (|8.24p in the proof of Theorem 16.11 below, we conclude that this proportion is 

/ dpr = bj, 

where bj is the period in (|2.12p . Thus, we easily deduce the following: 

Corollary 6.2. //, for a particular j = 1, . . . , N , Tj is singular, then there is a exists a subsequence 
{P n }neM, such that P n = QkRk, deg(Q fc ) = n k , where 

— VQ k — *-> PfItj, as n — > oo, n £ Af (6.10) 

and 




As stated in (iv) of Theorem 16.11 if a is a weak* cluster point of the measures {vp n }^ = i then: 
(a) supper C G and (b) the balayage of a onto T equals the equilibrium distribution pr- The 
following corollary shows that the equilibrium distribution is also obtained if weak* convergence 
and balayage are applied in the opposite order. 

Corollary 6.3. Let Bal {yp n ) denote the measure obtained by balayage ofvp n \c onto T while keeping 
vp n \<C\G unchanged. Then 

Bal {vp n ) — —> pr as n-> oo. 

6.2. Case studies. In this subsection we make more explicit Theorem 16.11 and its corollaries, and 
we illustrate them by means of a number of representative cases and examples. 

Case I: Two singular boundaries. 

Here N = 2 and p(<fj) = 1, j = 1, 2, for any two conformal maps ipj : Gj — > B. By Corollary 16 .1\ 
(3 equals the equilibrium measure of G and there exists, for each island Gj, a subsequence of vp n 
which converges to pr in a neighborhood of Gj. However, we do not know whether there necessarily 
exists a common subsequence for the two islands. 
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Case II: One singular boundary and one analytic boundary. 

Assume that Ti is singular and T 2 is analytic. Then in terms of two specific conformal maps 
ifj : Gj — ► B, j = 1,2: (a) ip\ has no analytic continuation beyond Ti, (b) (p 2 extends analytically 
as a univalent function to some domain containing G 2 . Since T2 is an analytic Jordan curve, it 
possesses a Schwarz function, which is given by 



S 2 {z) = ^\l/^ 2 (z)). 

In order to formulate a particular statement we assume further that tp 2 remains analytic and 
univalent throughout G 2 ,R'- This implies that gn(-,oo) extends by Schwarz reflection up to the 
level line L 2 ; see ()3.ip and the terminology in Example 13.11 Moreover, the domain 



' R' 



D2 '■= Ql,B! \ Q- 



R' 



is connected and is a domain of involution of the Schwarz reflection z 1— ► S 2 (z). 
Set 

E = G x u g 2 1 . 

It follows that the multi-valued function 

f $(z) ifz€C\G, 
$(z) := < . 7 — — \ v (6.11) 

1 ; \ l/&\S 2 (z)) if zeG 2 \g 2 ^. K ! 

is (locally) analytic in C \ E and (locally) continuous on C \ E. It also follows from the expression 
(|8,23p of p(K(-, z)) appearing in the proof of Theorem 16 .1\ by taking into account (|8.18p and (|8.19p . 
that 

1 if z € Gi, 

exp{ -g n {z,co)} if z 6 G 2 \ Q 2> ^ , ( 612 ) 



R' if z € Q 1 . 



The relations in (16. lip and (|6.12p yield at once, in view of Proposition 13. II and Corollary 13 .1\ the 
n-th root asymptotic behavior of {P n }^ =1 in C: 

( 1 if z e Gi, 
limsuplPn^)! 1 /™ = \ Mz)\ ifz,£C\E, (6.13) 

In addition, these relations provide more detailed information for the potential U@ of the canonical 
measure /3, and thus for the counting measures {^Pnj^i- 

Corollary 6.4. Under the assumption and notations of Case II, we have: 

' if zeG x , 

U l3 {z) = \ iog^ry -m{z,oo) ifzeC\E, ( 6 .14) 

In particular, 

(i) supp/3 = dE. 
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(ii) For any weak* cluster point a of {vp n }, supper c E, and 

U a (z) = U l3 (z), zeC\E. 

(iii) There is a subsequence N C N such that, with V any neighborhood of G\ or Q 9 1 not 
meeting the other island, 

vp n \v — * P\v, n —> oo, hGJV. (6.15) 
Hence, every point of dE = T\ U L 2 1 belongs to supper, for some weak* cluster point a of 
{"pJ%Li- 

The corollary is illustrated in the following example. 

Example 6.1. Bergman polynomials for G = G1UG2, with Gi the canonical pentagon with vertices 
at the fifth roots of unity and Gi = {z : \z — 7/2| < 2/3}. 

The zeros of the associated Bergman polynomials P n , for n = 80, 90 and 100 are shown in 
Figure In the same figure we also depict the critical line L^/ and the curve L 9 1 . Note that 

L 2 1 is simply the inverse image of with respect to the circle {z : \z — 7/2| = 2/3}. 




Figure 8. Zeros of Bergman polynomials P n of Example 16.11 for n = 80, 90 and 100. 



Case III: Two analytic boundary curves. This is the case N = 2, where both Ti and T2 are analytic 
curves. 

Example 6.2. Bergman polynomials for the union of the disks: G\ = {z : \z + 2| < 1} and 
G 2 := {z : \z-3\ < 2/3}. 

Let Si and 5*2 denote the Schwarz functions defined by Ti and T2- (Note that the Schwarz 
function for the circle {z : \z — a\ = r} is simply S(z) = r 2 / (z — a) + a.) Clearly, the Green function 
go, extends by Schwarz reflection to the set 

(6.16) 



D = {Gi, R/ \g 1 1 )u(g 2tR ?\g 2 o, 

' R' ' R> 



and the multi-valued function 
$(z) := 



$(z) ifzeC\G, 



l/*(Sj(z)) if zeGjXG^, i = 1,2, 



(6.17) 
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is (locally) analytic in C \ E and (locally) continuous on C \ E, where now 

e = g, i u g 2 i . 

As in Case II, the extensions of gci(z, oo) and &{z) lead to the expressions 

pW,,)) = { eXp{ -^'°° )} l* z % G E ) E > (6-18) 

hmsuplP^,)! 1 /™ = I 1*^)1 * Z ll^ (6.19) 

n—>oo ^ "^7 II ^ t fi, 

and in parallel with Corollary 16.41 to the conclusion 

„ | log — 3_ — Qn(z, oo) if z G C \ £7, 
[^( Z ) = J &ca PP yUV ; N (6.20) 

supp/3 = and that every point of dE = L l j_ U L 2 j_ attracts zeros of the sequence {-P n }^=i- 

Furthermore, since the unbounded domains C\E and 0, i coincide, it follows from (|2.14p . (|2.15p 
w 

and (|6.20p that the same is true for the potentials U@ and JJ^ 3E in C. Hence, the canonical measure 
j3 is the equilibrium measure of dE. Therefore, by applying Corollary 16.11 (ii) (with E in the place 
of G), we conclude that for j = 1,2, there is a subsequence J\f = Mj C N such that, with V any 
neighborhood of Q- i not meeting the other island, 

u Pn\v — -> ^8e\v, n^oo,n£j\f. (6-21) 




Figure 9. Zeros of Bergman polynomials P n of Example 16.21 for n = 140, 150 and 160. 

The zeros of the associated Bergman polynomials P n , for n = 140, 150 and 160 are shown in 
Figure [9l In the same figure we also depict the critical line Lri and the curves i and L 9 i . 

Note that L ■ i is the inverse image of La w with respect to the circle r,-, j = 1, 2. 

J'W ■" J 

Example 6.3. Bergman polynomials for the union of an ellipse and a disk. 

In Figure [10] we plot the zeros of the Bergman polynomials P n , for n = 80, 90 and 100 of an 
ellipse (domain Gi) and a disk (domain G^)-, in relative positions chosen to illustrate further the 
theory given above. To this end, let S\ and S2 denote the Schwarz function associated with the 
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(iii) 

Figure 10. Zeros of Bergman polynomials P n of Example 16.31 for n = 80, 90 and 100. 



ellipse, respectively the circle. The three ellipses pictured in Figure [TOl have all focal segment [-1,1] 
and canonical equation 



x2 y2 

a 2 + b 2 ' 



with a = 5/3, 6 = 4/3, in (i) and a = 5/4, b = 3/4, in both (ii) and (iii). 
For such ellipses the associated Schwarz function is given by 



Si(z) = (2a 2 - l)z - lab^z 2 - 1, 

and the focal segment [—1, 1] is reflected to the confocal ellipse x 2 / A 2 +y 2 / B 2 = 1, where A = 2a 2 — 1 
and B = 2ab. We denote by 

Di = {(x,y): ^ + ^<l}\[-l,l], 

the maximal domain of involution for the Schwarz reflection and by 7 the outer boundary of D±, 
i.e., 

x 2 y 2 
7 = {(*.?/) : a 2 + B 2=1} - 



Also, if G2 is a disk centered at z = zq, the reflection z 1— > S2 {z) is an involution on the domain 
D 2 = C\{z }. 

The situations illustrated in Figure [T0l represent the three possible relative positions between the 
loop Li^ri of the singular level set and 7: 

• Figure [10] (i) corresponds to the case that Li^> is interior to 7, 

• Figure [10] (ii) corresponds to the case that Li ri intersects 7, 

• Figure [10] (iii) corresponds to the case that the ellipse 7 is interior to L\^i. 
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: 



Figure 11. Zeros of Bergman polynomials P n of Example 16.41 for n = 80, 90 and 100. 

By specializing Theorem 16. II to this example, we can conclude the following: 

Case (i) is completely analogous to Example I6.2L That is, supp/3 = dE = L, i UL i and 

'W 'W 

every point of dE attracts zeros of the sequence {P n }^ =1 . More precisely, (3 = \xqe and for any 
j = 1, 2, there exists a subsequence M = M-j C N such that, with V any neighborhood of Q- i not 
meeting the other island, 

u P n \v HbeW, n -> oo, n £ J\f. 

In case (ii), the support of the canonical measure (3 consists of three parts: the inverse image 
Lr, i of L2 r' with respect to the circle T2, the reflection of L\ r> n D\ with respect to the ellipse 

T\ and the part [s, 1] of the focal segment [—1, 1] of the ellipse that lies exterior to this reflection. 
In addition, every point of supp/3 attracts zeros of the sequence {P n }%Li- 

Finally in case (hi), supp/3 = [—1, llUL, 1 . Thus C \ supp/3 has exactly two components and 

it follows from (vii) of Theorem 16.11 that there exists is a subsequence M C N such that 

v Pn ^ [3, n -> 00, n G N. (6.22) 

Case IV: One piecewise analytic non-singular boundary and one analytic boundary curve. 

Assume that T2 is analytic and T± is piecewise analytic and non-singular. By the latter we 
mean that any conformal map (p± : G\ — ► D has an analytic continuation to a neighborhood of G±, 
but this continuation is not univalent in any neighborhood of G\. This occurs, for example, if Y\ 
consists of circular arcs and/or straight lines and all its interior corners are of the form 7r/m, m > 2 
an integer. 

Example 6.4. Bergman polynomials for the union of the half-disk G\ = {z : \z\ < 1, He(z) > 0} 
and the disk G2 = {z : \z — 3| < 2/3}. 

In Figure [TTJ we plot the zeros of the Bergman polynomials P n of G, for n = 80, 90 and 100. In 
addition we depict: 

• The critical level line Lr> of the Green function (70(2,00). 

• The part of the reflection (we denote it by F^) of L\,ri with respect to Ti which lies in G\. 

• The inverse image L 2 1 of £2,/?' with respect to the circle T2- 
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Figure 12. Zeros of Bergman polynomials P n of Example 16.51 for n = 80, 90 and 100. 

By considering the symmetric and inverse images of the interior points of Gi with respect to 
the two arcs forming T\, in conjunction with the harmonic extension of the Green function inside 
G\ defined by the Schwarz functions of these arcs, it is not difficult to see that the support of 
the canonical measure (3 consists of three parts: the loop and two (symmetric) arcs that join 
together each one of the points i and — i with the nearest corner of T^. In addition, every point of 
supp/3 attracts zeros of the sequence {P n }^ =1 . 

Example 6.5. Bergman polynomials for the union of the symmetric lens domain G\ formed by two 
circular arcs meeting at —i and i with interior angles ir/4 and the disk G2 = {z : \z — 5/2| < 2/3}. 

In Figure IT21 we plot the zeros of the Bergman polynomials P n of G, for n = 80, 90 and 100. In 
addition we depict: 

• The critical level line Lri of the Green function gn(z, 00). 

• The part of the reflection (we denote it by T^) of Li r> with respect to T\ which lies in G\. 

• The inverse image L 9 1 of L2 r> with respect to the circle T2- 

As it is expected, identical conclusions to those of Example 16.41 regarding the properties of the 
support of the canonical measure (3 hold here. 

Case V: Three analytic boundaries. 

Example 6.6. Bergman polynomials for the union of the three disks G\ = {z : \z + 1| < 1/2}, 
G 2 = {z : \z - 2| < 1} and G 3 = {z : \z - 2i\ < 1/2}. 

In this example we have two critical Green level lines, Lr> and L_r», where R 1 = R2 = R3 and 
R" = R\. (See Figure CD which depicts the present example.) On setting 

E' = Q 1 U Go 1 and E" = Gi 1 , 
2 <W S 'W ~~rF 

we have 

( \$(z)\ if ze C\(E'UE"), 
limsup|P n (z)| 1 / n = <^ 1 iizeE', (6.23) 
( ir iizGE", 

where $(z) is the multi-valued function defined as in (|6.17j) . with j = 1, 2, 3. From (|6.23p and (|3.9p 
conclusions can be drawn about the canonical measure f3. In particular we note that supp/3 = 
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Figure 13. Zeros of Bergman polynomials P n of Example 16.61 for n = 80, 90 and 100. 
dE' U dE" = L 1 j_ U L 2 _i_ U L 3 j_ and that every point of dE' U dE" attracts zeros of the sequence 

{Pn}^Ll- 

In Figure IT31 we plot the zeros of the Bergman polynomials P n of G, for n = 80, 90 and 100. In 
order to illustrate the above observations regarding the zero distribution we also depict the inverse 
image L ■ j_ of La r. with respect to the circle Tj, j = 1, 2, 3. 

We end this section by noting that the critical level curves of the Green function depicted in 
the plots above, were computed by a simple modification of the MATLAB code manydisks .m of 
Trefethen [34] . The original code manydisks. m is designed for archipelagoes formed by circles; see 
also Remark 12.11 

7. An example: lemniscate islands 

Let G := {z : \z m — 1| < r m }, m > 2 an integer and < r < 1. Then G consists of m islands 
Gi,G 2 , ■ ■ ■ ,G m , where 

Gj contains e 2nji/m , j = 1, 2, . . . , m. (7.1) 

Let P n {z) = \ n z n +- ■ ■ denote the (orthonormal) Bergman polynomial of degree n for the archipelago 
G, and write 

n = km + s, < s < m — 1. 

By the rotational symmetry of G and the uniqueness of the Bergman polynomials it is easy to see 
that 

Pkm+s(z) = z s Q k)S (z m ), degQ k)S = k. (7.2) 

Then 

Pkm +S {z) := Pkm+s(z) = z s q Ks {z m ) = z km + s + • • • , (7.3) 

"■km+s 

are the monic Bergman polynomials. Our first result concerns the asymptotic behavior of the 
leading coefficient A n . 

Proposition 7.1. For each s = 0, 1, . . . , m — 1 there holds 

lim X km+S r kM \l—f^- = (7.4) 



BERGMAN POLYNOMIALS ON AN ARCHIPELAGO 29 

Remark 7.1. Note that r = cap(G) = cap(r), where as above T = dG. Thus the sequence 

A n cap(rr +1 ,/^^, ra £ N, 
V n + 1 

has exactly m limit points, -±rr, -^=j, . . . , ±, 1. 

In Table [77L1 we illustrate Proposition 17. II for the lemniscate depicted in Figure [T4l where m = 3 
and r = 0.9. More precisely, Table I7TT1 contains the computed values of the leading coefficients X n 
correct to 6 decimal figures, for n = 38, . . . , 52, together with the computed values of A„,r n+1 y^^py. 

As predicted by the theory, the values of A n r n+1 ^J^fi alternate, as n increases, towards to the 
three limits 

1/r 2 = 1.234567..., 1/r = 1.111111 ... , 1. 
The coincidence for the values of n = 38, 41, . . . , 50 is explained in the proof of Proposition 17.11 



n 


An 


\ r n+l 1 7r 


38 


214.535664 


1.000000 


39 


305.078943 


1.263740 


40 


305.314216 


1.124276 


41 


305.396681 


1.000000 


42 


433.231373 


1.261795 


43 


433.526043 


1.123400 


44 


433.629077 


1.000000 


45 


613.834469 


1.260094 


46 


614.205506 


1.122633 


47 


614.334958 


1.000000 


48 


868.011830 


1.258593 


49 


868.481244 


1.121956 


50 


868.644692 


1.000000 


51 


1225.297855 


1.257261 


52 


1225.894247 


1.121355 



Table 7.1. Illustrating Proposition 17.11 for the lemniscate case m = 3 and r = 0.9, 
for n = 38, . . . , 52. 



Proposition 7.2. The following representations hold for the monic polynomials Pkm+s( z ) : 

p km+m . 1 (z)=z m - 1 (z m -l) k (7.5) 
and for s = 0, 1, . . . , m — 2, we have for k sufficiently large, 

= *k + i,sW ~ M (7-6) 
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where w = (z m — l)/r m and ir n>s (w) is the monic polynomial of degree n in w that is orthogonal on 
the circle \ w\ = 1 with respect to the weight 

\dw\ 



l w + II 2 " 



_2__ 2s 

m m 



(7.7) 



Remark 7.2. The representation formulas (|7.5p and (|7.6p have the same form as those found by 
Miha-Dfaz p3], who studied the simpler case when r > 1, i.e. when G consists of a single island. 



In our proof we utilize the following lemma that relates 'weighted' Bergman polynomials on the 
unit disk to Szegd polynomials on the unit circle. This result is somewhat implicitly contained in 

mi. 



Lemma 7.1. Let t n (w) = w n + • • • be the monic polynomial orthogonal with respect to the weight 
\dw\/\~fw + 1| T on \w\ = 1, where r is real, r ^ 2,4,... ,2n, and \~f\ < 1. Let (3 n (w) = w n + ■ ■ ■ 
be the monic polynomial orthogonal with respect to the weight dA(w)/\'yw + l| r over the unit disk 
\w\ < 1. Ift n (—j) / 0, then 

(w + j)P n (w) = t n+l (w) - ^Z^ t n (w). (7.8) 

Our next result describes the fine asymptotics for the monic Bergman polynomials. 
Proposition 7.3. Let 

T :=2 - — -—, s = 0, l,...,m- 1. (7.9) 
m m 

Then for \z m — 1| > r 2m ,z m - 1 / -r 2m , the monic Bergman polynomials satisfy for each s = 
0, 1, . . . , m — 1 

iim (^i±^Y'\ (710) 

fc-^oo Z s (z m - l) fc V. Z m - 1 / 



where the branch of the power function on the right-hand side of {1. 10 ) is taken to equal one at 
infinity, and the convergence is uniform on compact subsets. 

Furthermore, for each j = 1, 2, . . . ,m and z S Gj with \z m — 1| < r 2m , we have 

(-l) k+1 k 2+T / 2 e 27ry( s +i)/m z m-l rr ( r / 2 ) s in(r7r/2) 

/or eac/i s = 0, 1, . . . , m — 2, the convergence being uniform on closed subsets. 

Observe that the lemniscate \z m — 1| = r 2m is the reflection of the lemniscate \z m — 1| = 1 in the 
bounding lemniscate of G. 

Remark 7.3. From the first part of Proposition 7.3 we see that the Bergman polynomials for G 
have no limit point of zeros in \z m — 1| > r 2m other than at z = 0. Furthermore, from the second 
part of the proposition, we deduce that, except for the subsequence (17. 5p . there are no limit points 
of the zeros of P n {z) in \z m — 1| < r 2m . Consequently, the only limit points of zeros of such P n (z) 
are at z = or on the lemniscate \z m — 1| = r 2m . 

In Figure [TU we plot the zeros of the Bergman polynomials P n , for n = 50, 51 and 52, of G := 
{z : \z 3 — 1| < 0.9 3 }. In each plot, we depict also the defining lemniscate T = {z : \z 3 — 1| = 0.9 3 }, 
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Figure 14. Zeros of the Bergman polynomials P n for the lemniscate case m = 3 
and r = 0.9, for n = 50, 51 and 52. 

the reflection {z : \z s — 1| = 0.9 6 } of {z : \z 3 — 1| = 1} in T and, for the cases n = 51,52, 

z 3 _ i _|_ q g6 

branch cuts for the Schwarz function S(z) = ( 5 ) 1//3 of F. 

z 6 — 1 

As a consequence of Proposition 7.3 we have the following: 
Corollary 7.1. There are precisely two limit measures for the sequence {^p n }^ = i; namely 

^ m 

~ ^ S *i ' z i = exp(27rij/m), 
i=i 

and the equilibrium measure for the lemniscate \z m — 1| = r 2m , which is given by the formula 

l^lm— 1 
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8. Proofs 

The present section is devoted to the proofs of the results stated earlier in the article. 

Proof of Lemma 13. 1L That T is analytic is clear, since u is real analytic and Vu ^ 0. 

All of D is filled with integral curves of the gradient Vu. These are disjoint and have no end 
points in D since Vu 7^ 0. Hence they all end up on dD (an integral curve cannot be closed since 
u is single- valued and increases along it). These integral curves are at the same time level lines of 
any locally defined harmonic conjugate of u. 

Given z G D we want to define the reflected point S(z) using only u. Assume for example that 
u(z) < 0. By the maximum principle, \u\ < c in D, so actually — c < u(z) < 0. There is a unique 
integral curve 7 of Vw passing through z, and u increases along 7 with limiting value +c as 7 
approaches dD. Thus there is a unique point w G 7 at which u{w) = —u(z). In terms of this we 
define 



S(z) = w. 

The above procedure defines a function S(z) in D. To see that S(z) is analytic, note that, in 
some neighborhood of 7, u has a single-valued harmonic conjugate u* and that 7 is a level line of 
it*. The function / = u + iu* is analytic in a neighborhood of 7, with /' 7^ 0; hence / can be used 
as a new complex coordinate near 7, or u and u* are new real coordinates. In terms of these, the 
reflection * ^ S(z) just defined is given by 

u + iu* 1— > — u + ut* 



or /(z) i — > —f(z). This gives 



S(z) = f-i(-f(z)), 

which proves that is analytic. It is also immediate that S(z) = z on T, so that S is indeed a 
Sch war z function of T. □ 
Proof of Lemma 13.21 According to Theorem 3.2.3 of [27], one criterion for gL4|g to belong to 

the class Reg is that 

lim WPJ 1 ^ = 1; (8.1) 

n— >oo Cr 

note that f2 is regular with respect to the Dirichlet problem [20^ p. 92]. (Here and in the sequel 
|| • || means the sup norm on the subscripted set.) 

The argument given in the proof of Lemma 4.3 of [18], when separately applied to each of the 
Jordan regions Gj yields 

limsup||P n ||' /n < 1, j = l,2,...,iV. 

ill /n 11I In 

Consequently, limsupj^oo ||P n ||^ < 1- But liminf^^oo ||P n ||^ > 1> since ||Pn.||L2( G ) = 1 for all 
n, and so (|8.ip follows. □ 

8.1. The extremal problems. We use V n to denote the space of complex polynomials of degree 
n. Recall that K n (z, £) denotes the n-th finite section of K(z, £) 



k=0 
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K%{z,Q :=^P fcii (C)P fc j(z), 



and similarly set 

k=0 

where 

P n ,j{z) = X njz n + ■■■ , X n>j >0, n = 0, 1, 2, . . . , 
are the sequences of the Bergman polynomials associated with Gj, j = 1, 2, . . . 
Lemma 8.1. For any (gC, 



max 



VX(C.C), n = 0,l,.... 



pdVn \\P\\ L 2 {G) 

Proof. Since for any p G V n and Q G C 

p(C) = (p,^»(-,C)>, 

it follows 

|P(C)I < IHI^G) ll^n(-,C)|k 2 (G) = l|p||L 2 (G)V^n(C,C)- 

Hence 

INU 2 (G) 

with equality if p(z) = cK n (z, (), for some constant c / 0. 
Obviously 

WpWmGj) < IIpIIl 2 (g), j = 1,2, . . . ,iv, 

therefore for n = 0, 1, . . ., 

maoc^L^max^L, ; = 1, 2, . . . , N, 



or 



K^((,0>K n ((,0, j = l,2,...,N, (SC. 



Furthermore, since for any £ G Gj, 

K% (C, C) < K G > (C, C) = tf(C, 0, 3 = 1, 2, . . . , JV, 
it follows from (|8.2[) that 

1 1 1 

> ; > t=^==f , j = l,2,...,iV. 



v^U^C) yjK°>(c,c) v^CCy 
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8.2. Proof of Theorem 14. 1L The estimates from above require only a C 2+Q -smooth boundary 
and are based on comparison with corresponding estimates for the arc-length measure \dz\ and the 
Szego orthogonal polynomials. To this purpose, we compare the two extremal problems 

m 2 n {G,dA):= min / \z n + a n _iz n_1 • • • + a \ 2 dA(z), n = 0,1,2,..., (8.4) 

ao,...,On_i Jq 

and 

,2/"n „l-7~h . / I jn i „ jn,— 1 



m z n {T,p\dz\) := min / \z n + a n -iz n 1 h a \p(z)\dz\, n = 0,1,2,..., (8.5) 

where p is a positive smooth function on V. Recall from (12. ID that 



m^G,dA) = -2= / (8.6) 

A n JG A n 

where 

P n (z) = X n z n + --- , A n >0, n = 0,l,2,..., 

are the Bergman polynomials of G. 

The asymptotic properties of m n (T,p\dz\) have been established by Widom in [38^ Thm 9.1]. In 
particular, the next estimate for p = 1 and some constant C > follows from Theorems 9.1 and 
9.2 of [38]: 

m 2 n (T,\dz\) >Ccap(r) 2 ". (8.7) 
On the other hand, Suetin's lemma (Lemma 13.51 above) applied to each island separately gives 

m n (G,dAf= [ \?MfdA>-^- [ \?M\*\ dz \>-?-rn n {TMz\)\ 
Jg *n n + 1 J T \ n n + l 

where C > is a another positive constant. 

Combining the above two estimates we conclude 

cap (r) re 



m n (G,dA) >C 



n 



which yields the upper inequality in Theorem 14.11 

For estimates from below we require analyticity of the boundary. The main technical aid is 
provided by a family of polynomials uj n constructed by Walsh in [35], which we thereby refer to as 
Walsh polynomials. 

Lemma 8.2. Assume that each Tj, j = 1,2, ... ,N, is analytic. Then, there exists a sequence of 
monic polynomials ui n {z) = z n + ■ ■ ■ , n = 1, 2, . . . , with all zeros on a fixed compact subset E C G, 
and a constant C such that 

IKIIl2(g) < 4= ca p( r ) n - (8-8) 

v n 

From this we deduce the lower inequality in Theorem 14. II 
Corollary 8.1. // each Tj, j = 1, 2, . . . , N , is analytic then 
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Proof of Lemma 18.21 Since each Tj, j = 1,2, . . . , N , is analytic, the Green function gn(z, oo) 
extends harmonically across dG by Schwarz reflection. Choose first a number < r < 1 such 
that - < R' (see Subsection 12.41 for the definition of R') and such that gn(z, oo) extends into each 
component of G, at least to the negative level log r. Since gn(z, oo) has no critical points in Qri \ G 
it follows that the extended Green function has no critical points in D = Qi \ Q T = Q T \ 

T 

The latter open set has N components, each of which is a domain of involution for the Schwarz 
reflection (see Lemma l3.ip . 

Now choose a number p in the interval 

r < p < 1. 

For any R> p, 

gn R {z, oo) := gn(z, oo) -log R, (8.10) 

is the Green function of £Ir with pole at infinity. Hence, 

cap(L R ) = iicap(r). (8.11) 

Choose the compact set E C G in the statement of the lemma to be E = L T . By a theorem of 
Walsh [35] (see also [2TJ p. 515]), there exists a sequence of monic polynomials uJ n (z) = z n + ■ ■ ■ , 
n = 1, 2, . . . , with zeros approximately equidistributed with respect to the conjugate function of 
gn(z, oo) and such that 

1 C 

\gn T (z, oo) +logcap(L r ) log |w„(z)|| < — in VL p . (8.12) 

n n 

Note that gn R (z, oo) + log cap (L^) is independent of R, hence in (|8.12[) r can be replaced by 
any number R > r. For z £ Lr and R> p this gives 

1 C 

|logcap(LR) log 1^(^)11 < — , 

n n 

or after exponentiating and using (|8.1ip 

e~ c < <e c (zeL R ,p<R<oo). (8.13) 

In particular, from the maximum principle, 

\uj n {z)\ < CR n cs^{T) n (z € G R , p < R < oo), (8.14) 

for another constant C. 

Next we estimate the L 2 (G)-norm of u n . On decomposing 

uo n \ 2 dA= I \uj n \ 2 dA+ / \oj n \ 2 dA, 

G JG P JG\G p 

the first term can be directly estimated by means of fj8. 14j) : 

/ \u n \ 2 dA < Cmax \to n (z)\ 2 < Cp 2n cap(T) 2n . 
JG P 2£L p 

For the second term we foliate G\G P by the level lines Lr of gn(z, oo), or (^(z)! = exp[gfj(z, oo)], 
and use the coarea formula. Since Vgn(z,oo), and hence V|<I > (z)|, is bounded away from zero on 
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G \ G p we obtain by using once more (|8.14|) 



\u n \ 2 dA = ^ f ^ri.Hi, W\dR 
G\G P ' J P Jl r |V|*(*)|| 1 1 

c 



[ max \uj n (z)\ 2 dR < Ccap(r) 2n / R 2n dR 

J O Z ^ L R Jo 



<- ^^^^ 

for various positive constants C. Thus altogether we have 



/ M 2 ^<C(p 2n + -)cap(r) 2r \ 

Jg n 



and since p < 1, this gives (18.81) . □ 
The corollary is an immediate consequence of the lemma and the definition of X n : 



-3- = m n (G, dA) < ||w n ||L 2 (G) < C 



8.3. Proof of Theorem 14.21 We turn now our attention to the problem of determining the rate 
of convergence of A n J as compared to A n . The solution will obviously depend on a set of numerical 
constants which reflect the global configuration of G. 

In the case of a single island TV = 1 we have A^ 1 = A n , hence both (|4.7p and (|4.8p hold trivially 
with m = 1. For the case N > 2, we assume that Tj is analytic, for some fixed j G {1, 2, . . . , N}. 
Let <Y denote the characteristic function of Gj in G and set 

. , II^pIU^g) 

7„ := mf — — (8.15) 

pePn \\p\\l2(G) 

(Note that 1 1 p 1 1 x, 2 = |b||i2 (Gj ), hence < 7„ < 1.) 

By considering the Bergman polynomial P n j of Gj, as a competing polynomial in (|8.15|) and using 
Carleman asymptotics (Theorem 13. ip for P n j in G\Gj in conjunction with the fact |$j(z)| > |$(z)|, 
z E CI (subordinate principle for the Green function; see e.g. |2U|, p. 108]), we conclude that there 
exist constants C > and R > Rj (> 1) such that, for any n G N, 

— > 1 + C v 7 ^ R n - 

In 

Hence for large values of n, 

In < a n , 

where < a < 1. Since X has an analytic continuation up to Lri in Q, it follows from Walsh's 
theorem of maximal convergence [36^ Thm IV. 5] that for any n G N, there exist a constant m > 1 
and a polynomial g m ( n ) G "P m („), where m(n) = mn, with the property, 

hm(n) ~ X \\g < 7n- (8-16) 

Then we have: 
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Lemma 8.3. Assume that Tj, j S {1,2,... , N}, is analytic. Then for any n G N 



Af J ((,()< ~r~2 \/K n+m{n) ((,(), C^Gj. 
Proof. Take £ E Gj and let h £ V n be an extremal polynomial for 

|p(C)l 



max ■ 



Then from Lemma 
It holds, 

because from fl8.16|> . 
Also 



peVn | |P 1 1 £,2 ( G ) 



\\Xh\\ L 2 {G) - 



\(Xh)(()\<-^—\(q m{n) h)(0\, 



(l- 7n )^(C) < |«m(„)(C)|- 



\\Qm(n)h\\L 2 (G) < \\Xh\\ L 2 (G) + || (A? - g m(n ) )h\\ L 2 ( G ) 

< ll^ll^(G)+7n||^|k2(G) < 2||^|| i2(G) , 

where in the last inequality we made use of the defining property of 7 n . Finally, 

\(Xh)(0\ < 2 |(g m(n )/i)(C)l 



I^^IIl 2 (G) 1 - 7n |km(n)^lli 2 (G) 

. 2 |/(C)| 
< max 



1 -In f£P n+m(n) ||/|| L 2 (G) 

and the result follows from Lemma l8.11 □ 
This yields Inequality (14. 8j) in Theorem [421 The other inequality (14. 7ft follows immediately from 
(El. 



8.4. Proof of Theorem 14.31 Keeping in mind Lemma 13.31 it is clear from its definition that 
the functions A n converge uniformly on compact subsets of G to A. By imposing analyticity of 
the boundary, we will be able to estimate jointly the rate of convergence of A n (z) on T and in a 
neighborhood of T in the interior. In view of the reduction to a single island established in the 
previous subsection, we will assume in the first part of the proof that N = 1. In order to simplify 
further the notation, we will simply write G = G\. $ = <I>i and so forth. 

Thus, we deal now with a Jordan domain G with analytic boundary T. The normalized external 
conformal mapping <I> analytically extends to the level set Q p , with p < 1. According to Theorem l3.ll 
the Bergman orthogonal polynomials satisfy: 



\l T ^${z) n &{z){l + A n {z)}, z£G\G p , 



K(z, z) = 1 V _ ;| M + R{z). 
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where A n (z) = 0((^) n ), whenever z £ T r , and p < r < 1. Fix a z £ G \ Q p and denote t = \§(z)\ 2 . 
Then 

K n (z,z) = J2\P k (z)\ 2 = 1 -^- J2( k + l)t k + R n (z) (8.17) 

k=0 k=0 

\&(z)\ 2 l-{n + 2)t n+l + ( n + l)t"+ 2 

Similarly, 

|$ '(X)| 2 i_ 

The convergence of R n {z) to R(z), for p 2 < r 2 < i < 1, is uniformly dominated by a convergent 
geometric series. 

In view of (|4.4j) we set A(z) = for all zeT, Since 

< A n (z - A(z) = 1 1 

y/K n (z,z) yJK(z,z) 

we are led to the estimate 

AJz) - A(z) < C(l - t)[ ; 1 - 11. 

" ' VI - (" + 2)i n+1 + (n + l)i"+ 2 

In its turn, elementary algebra yields: 

, _ . r 1 , 1 (n + 2)t" +1 - {n + l)t n + 2 



■y/l-(n + 2)* n + 1 + (n + l)i n + 2 \/£fc=o( fc + 1 )* n 1 + V 1 = ( n + 2 )* n+1 + (™ + i)^ 2 

- W 2 A /l + 2 + ... + (n+l) n + 1 J 

< Ct n / 2 (1 -t+ -), 
n 

which implies Inequality (|4.9p in Theorem 14.31 since for z near T: 

i - |$(z)| 2 x i - 1$(«)| x dist(>,r). 

Using (|8.17p . which holds for z £ F with R n (z) = 0(n 2 y/np n ), we derive easily (|4.10p . which is 
the limit of the exact form of (14.91). 



We resume now our general assumption G = Uj =1 Gj and we turn our attention to deriving 
(|4.1ip . The lower bound emerges at once by combining (|4.10p with (|4.7p . To obtain the upper 
bound we apply (|4.8p to A/%(z), for large k, with k = [k/m]m + r, where < r < m — 1, and [fc/m] 
is the integral part of the fraction, and then we use again (|4.10p . 

In order to estimate A n in the exterior of G we employ the Walsh polynomials: From Lemma [8 -H 

IHIl 2 (g) 



A, 



mm 



P &v n \p{z)\ 
and therefore, 



\u n (z)\ - 

where we made use of Lemma 18.21 and f|8. 13|) . 
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Finally, the lower estimate for A n (z) for z exterior to G is directly derived from the upper 
estimates for the orthogonal polynomials appearing in Theorem I4.4L □ 

8.5. Proof of Theorem 14.41 Our aim is to derive estimates for P n (z), for z in the exterior of 
the archipelago. To do so, we assume that every curve constituting T is analytic and we rely, once 
more, to the Walsh polynomials uj n . 

P (z) 

We fix a positive integer n and consider the rational function ^ > whose poles lie in a compact 
subset of G and which vanishes at infinity. With z ^ G, Cauchy's formula yields: 

P n (z) -1 f P n (()d( 



uj n+ i (z) 2m J r ui n+ i (C) (C - z) ' 

whence, from f|8.13j> . 

IPMK C lip II 

l ^ n{Z)l ~ dist(z,T) cap(r)"+i l|jPre|lil (r)> 

where || • ||l1(t) denotes the L -norm on T with respect to \dz\. 

Since the L -norm is dominated by a constant times the L 2 -norm, Lemma [3.51 gives H-PnHi^r) — 
C y/n and one more application of (|8.13p yields 

(In the above we use C to denote positive constants, not necessarily the same in all instances.) 

In order to obtain the estimates from below, we have to restrict the point z to the complement 
of the convex hull Co(G). On that set, including the point at infinity, the sequence of rational 
functions R n (z) = Pn (^) ca P(r) ^ as no zeroS; anc [ by the above estimate, it is equicontinuous on 

compact subsets of U = C \ Co(G). Thus {R n }^L forms a normal family on U and the possible 
limit functions are either identically zero, or zero free. The normalization at infinity was chosen so 
that, in view of (|4.2p and (|8.13p . inf ne pj R n (oo) > 0. Thus, every limit point of the sequence R n is 
bounded away from zero, on compact subsets of U. 

8.6. Distribution of Zeros. 

Proof of Theorem 16.11 To prove (i) we need to figure out the general structure of p(K(-,z)). 
We have already remarked, cf. (|3.9j) . that for ( £ Gj, 

p(K(;C))=mm{R j ,p(K G ^;0)}. 
Recall (]2.17p . that is, in terms of any conformal mapping (fj : Gj — > B, 



7T [1 - <Pj(z)<Pj(0 

Conversely, if (given ( £ Gj) ipj is chosen so that </?j(C) = 0, then 

iPj( Z )=^= r K G i(t,c)dt. 

<fi(Q J c 
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Hence, for a general <pj, 

<Pj(*)-<Pj(0 _ ^(1 — lvj(C)l : 



j Z K G >{t,Qdt. 



i-^(^(C) ^(C) 

It follows therefore that, given a £ € Gj and a simply connected region .D with Gj C D C Gj,Rj, 
K j (z, C) has an analytic extension to D as a function of z if and only if tpj(z) has a meromorphic 
extension to D and does not attain the value l/<^j(£) there. 

We introduce a meromorphic version of the function p defined in (|3.7h by setting, for / mero- 
morphic in G, 

Pm(f) '■= sup{-R > 1 : / has a meromorphic continuation to Gr} . (8.18) 

Next we extend each ipj to all G by setting (fj = in G \ Gj . Clearly the so extended ipj cannot 
be meromorphic in Gj,R for any R > Rj, hence 

1 < Pm(fj) < Rj. (8-19) 

(This is vacuous statement if N = 1, thus we simply set R\ = +oo in such a case.) The largest R 
for which ipj does not take the value l/tpj(() in Gj,R is inf{|$(</?j|g ( . } (V^.? (C)))|} 1)) where 

the infinmum is taken over all points in the preimage tpjlg 1 { (l/tpj(Q), which is a subset of 

Gj, Prn (ipj) \ Gj. (We assign the value +oo for the infimum of the empty set.) 
Putting things together we get, in view of (|8. 19|) . 

p(K(.,C))=min\p rn (<p j ), inf{|$(rf .(1/^(0)) \}\ , ( € Gj, (8.20) 

or, by taking the logarithm, 

logp(A-(-,0) = min(logp m (^-). ™f{9n(<Pj\g 1 , (l/^(0),oo)}) , ( e G r (8.21) 

This may look messy, but in principle it means that we have expressed log p{K{-, £)) as the infimum 
of some harmonic functions. This is the basic argument telling that log p(K(-, £)) is superharmonic 
as a function of £ in Gj . 

Now, if (fj has a singularity on Tj, then p m ((pj) = 1 and p(K(-,()) = 1, £ € Gj. In the 
complementary case, i.e., if y?,- has an analytic continuation across Tj, then for any £ £ Gj, 
(fjlg 1 ( ,)(VVj(C)) i s either void or it defines a (possibly) multi-valued reflection map in Tj, 

i.e., the conjugate of a (possibly) multi-valued Schwarz function of Tj. By our assumption that 
the infimum of the empty set is +oo, we only need to concentrate on the latter case. Denoting 
fj\g- (1/Vj(0) by 5j >mu iti(C) we can write (I8.2ip somewhat more handily as 

,multi 

(C),oo)}}, C £ Gj, (8.22) 

where the infinmum is taken over all branches of £j )mu iti(C)- O ne step further, this reflection map 
gives a multi- valued analytic extension of the Walsh function $ into Gj: 



$multi(C) = V$ Si,multi(C) , C € Gj 
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(where we have used hat to emphasize the analytic extension) . Inserting the latter into (I8.20p gives 
the following, more direct, description of p(K(-,Q)): 



the infimum is taken, again, over all (local) branches. 

In order to make the above considerations more rigorous we take (|8.2ip as our starting point. 
We first treat the case N > 2, which is somewhat simpler because in this case (|8. 19j) gives an upper 
bound for logp(K(-, £)) in (|8.2ip . Let £ G Gj. Then l/cpj(Q is outside the closed unit disk, and 



is larger, than near Tj. This means that only finitely many of the points in the preimage will 
be serious candidates in the competition for the infimum in (|8.2ip . We may also vary £ within a 
small disk, compactly contained in Gj, and there will still be only finitely many branches of the 
multivaled analytic function involved, when forming the infimum. Within such a disk there 

will also be only finitely many branch points (where two or more preimages coincide). 

Thus, locally away from the mentioned branch points, log p{K{-, £)) is the infimum of finitely 
many harmonic functions, hence is continuous and superharmonic. At the branch points log p(K (•, Q)) 
is still continuous, and since the set of branch points is discrete (in Qj tPm Up) \ Gj) they make up a 
removable set for continuous superharmonic functions; see e.g. \20\ Thm 3.6.1]. It follows, therefore, 
that log p(K(-, £)) is superharmonic (and continuous) in all Gj. 

We apply now the above inferences to h(z) = — log p(K (■, z)) , for z G G. If Pmifj) = 1) f° r 
some j, then h(z) = 0, for z E Gj, hence the transition across Tj to gn(z, oo) is continuous and 
subharmonic. If Pmi^j) > 1 an d <Pj remains univalent in a neighborhood of Gj, then it is easy 
to see that h(z) defines the harmonic continuation of gu(z, oo) across Tj (in fact, Tj turns out to 
be analytic and thus SVmuiti is the associated ordinary single- valued Schwarz function). Finally, if 
p m {ipj) > 1 but (fj is not univalent in any neighborhood of Gj then locally, away from finitely many 
branch points on Tj, h is still the ordinary harmonic continuation of gn(z, oo). At the branch points 
h is still continuous and the set of branch points is too small to affect the overall subharmonicity. 
Hence, in all possible situations h(z) = — log p(K (■, z)) is continuous and subharmonic in G. 

Therefore, we have established so far that in the case N > 2, h is subharmonic (and continuous) 
in C and since it coincides with the Green function in f2, (3 is a positive measure, with support 
contained in G. Moreover, from Gauss' theorem (see e.g. [24^ p. 83]), and the singularity of the 
Green function at infinity, we have for any R > 1: 



Hence /3 is a unit measure and this completes the proof of (i), for N > 2. 

In order to derive (ii), we observe that the Riesz decomposition theorem for subharmonic func- 
tions applied to h in C (see e.g. [201 P- 76]) gives, 



where v is harmonic in C. Then, by considering the expansions near infinity of U^(z) and h(z) = 
ga(z, oo), we see that v(z) = —log cap (r), which yields (16. 3p . Relation (16. 4p is an immediate 




(8.23) 





(8.24) 



h(z) = -U p {z)+v{z), z G C, 
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consequence of ()6,3p the fact that h coincides with the Green function in 0, in conjunction with 
the relations ([23jl - ([2?7|l . 

When N > 2, U@ is bounded from above because of (|8. 19j) : 

< log 3 \^- < oo. 

cap (1 j 

Statement (iii) of the theorem is just a juxtaposition of Proposition 13.11 and Corollary 13.11 along 
with (jO|) . 

As for (iv), C is nonempty by general compactness principles for measures and the known fact 
that all counting measures up n have support within a fixed compact set; see Remark l3.1[ Let a G C. 
Then there is a subsequence M = N a C N such that 

vp n — — > a, n — ► oo, n G A/". (8.25) 

Using the lower envelope theorem [24^ Thms 1.6.9 ] and (]6.6p we get 

CT(z) = liminf U Vp -(z) > lim inf U Vp " (z) = U p (z), (8.26) 

n— >oo n— >oo 

where the first equality holds only quasi everywhere in C. However the relation between U a and 
U@ persists everywhere in C, since both members are potentials. 

Let T> be any component of C \ supp/9. Applying the minimum principle to u = U u — U@ > 0, 
which is superharmonic in T>, gives that either u > in all T> or u = in all T>. Since u vanishes 
at oo (recall that a and (3 are unit measures) it follows that it vanishes in the entire unbounded 
component of C \ supp/3. From this and the observations above follow all parts of (iv). 

Turning to (v), let 

U = lsc(inf IT 7 ). 

By (iv), U 13 <U in C. To prove the opposite inequality, choose an arbitrary point z G C. Then 
there is subsequence M z C N, such that the liminf in (16. 6j) is realized at z, i.e. 

lim U up ™(z) = U p {z). (8.27) 

n— >oo 

By weak* compactness there exists a further subsequence M' z C M z and a measure a = a z G C such 
that 

vp n — —> a, n — > oo, n G M' z . (8.28) 
Then, by the principle of descent (see |24} Thin 1.6.8]) and (|8,27p . 

U p {z) = liminf U""»(z) > U a {z). (8.29) 

n— >oo 

Since z G C was arbitrary, 

C/' 3 > inf J7 CT in C, 
o-ec 

by which > C/ follows in all C. 

To finish the proof of (v), we let again D be a component of C \ supp (3. By choosing above z £ T> 
we get a measure a = a z G C with U a {z) = U^(z) (since equality necessarily holds in (|8.29p ). Thus 
U a = U 13 in all V because, as we have already proved, the other alternative would be U a > U@ in 
all V. 
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Regarding (vi), if C consists of only one point, say a, then U 13 = U a by (v), and from the unicity 
theorem for logarithmic potentials (see |24} Thm II.2.1]) we must have j3 = a. Clearly, the full 
sequence must converge to 0, because otherwise one could extract a subsequence converging to 
something else, which would be a different element in C. 

The assertions in (vii) are easy consequences of (iv) and (v): Since, for any a G C, U a = U" in 
the unbounded component of C \ supp/3 we get in the case of (a) plus (b) that (for any a G C) 
U a = U@, almost everywhere with respect to the area measure in C. This and the unicity theorem 
yield = a G C. In the case of (a) plus (c), there exists (by (v)) at least one o G C satisfying 
U a = XjP in the bounded component of C \ supp 0, and for this a we have the same conclusion: 
U a = U@ almost everywhere in C and, as above, = a G C. 

So far we have assumed that N > 2. Let us indicate the modifications needed for N = 1. 
Equation (I8,21h may be written 

logp(K(-X))= hm min 

M— > + 00 *3,PmWj) J 

that is, by introducing an auxiliary upper bound M, which finally tends to infinity. Before passing 
to the limit we can work with the corresponding quantities 

h M = sup{h, -M}, M = ^-A/iju 
(etc.) as before. Since a decreasing sequence of subharmonic functions is subharmonic, h = lim Km 

M—*oc 

will be again subharmonic. It is however not clear that it will be continuous, only upper semicon- 
tinuity is automatic. If p m ((pj) < oo, then the bound M is not needed and everything will be as in 
the case N > 2. So assume p m (<Pj) = °°- This means that ipj is meromorphic in the entire complex 
plane and hence (|8.2ip reads 



logp(tf(.,C)) =inf{ ff n(^|c 1 (lM(0),oo)}, C G G 5 . (8.31) 

Problems concerning the lower boundedness and continuity of h could conceivably occur at points 
C G G at which the inverse image above is either empty or is an infinite set. The first case can, by 
Picard's theorem, occur for at most two values of Q G G. At such points the infimum in (|8.3ip is 
+oo, and hence h(£) = — oo. In particular, h will not be bounded from below, but it will still be 
subharmonic and upper semicontinuous. Moreover, it will be continuous at all other points, which 
is enough for the reasoning in the proof (above) of (iv), where we used the continuity of h (or U@). 

The second conceivable problem, that Lpi\^{\/ Lpi{C)) is an infinite set, presents no actual dif- 
ficulty because the only cluster points can be at infinity, hence all but finitely many branches of 
99i|^ 1 (l/(/9i(C)) will be ruled out when taking the infimum in (I8,31j) . □ 

Proof of Corollary 16.11 As already remarked, the boundary curve r, is singular if and only if 
p m (ifj) = 1, which by the proof of the theorem (e.g., Equation (|8.21jl ) occurs if and only if h = 
in Gj. This, in view of (|6.3p . is equivalent to 



Also from (16.3 
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It follows that is harmonic in Gj,R \ Tj , thus supp [3 C Tj . It also follows that the logarithmic 
potentials of (3 and fir coincide in the domain Gj,Rj, hence the equation @\q. = Uric holds as a 
result of the unicity theorem (see e.g. [Ml P- 97]). This proves the equivalence of (i) and (ii). 

By assertion (v) of the theorem, there exists a a G C such that U a = U@ in (= T>). The 
equation persists on Tj, because of the continuity of logarithmic potentials in the fine topology and 
in view of (|6.7p . it also holds in any neighborhood of Gj not meeting the other islands. Thus, from 
the unicity theorem a = f3, in such a neighborhood. As a is a cluster point of {^p n }, we conclude 
that (iii) follows from (ii). 

If (iii) holds then by selecting a further subsequence we conclude a\y = nr\v, f° r some a G C. 
Then U a = U^ r in V, which in view of (J6T1]) and (pr?) yields the relation = U^ F in V. Therefore 

Proof of Corollary D Set fi n = Bal (u Pn ). Then 

supp fin C C \ G, (8.32) 



= in 0. (8.33) 

Let fi be any weak* cluster point of {fin} and let M C N be a subsequence with \i n — n G M. 
By refining we may assume also that up n — —> a, n G J\f, for some measure a. Then in view of 
dSMD we have = in 0. 

On the other hand, U a = in VL by Theorem 16. 11 thus = in VL. But [/ Mr is harmonic 
in Q, \ {oo} and supp /U C C \ G by (|8.32|) . hence supp (iCT, Now Carleson's unicity theorem |24[ 
p. 123], shows that \i = fir- Since fi was an arbitrary cluster point of fi n it follows that fi n — — > /.tr 
for the full sequence. □ 

Proof of Corollary 16.41 The expression for U 13 follows immediately after uploading (I6.12p into 
Theorem l6.1l (ii). From this expression and the unicity theorem for logarithmic potentials we gather 
that supp/3 must be contained in dE. To show that eventually supp/? = dE we can argue as in 
[TrJl pp. 215-216]. That is, by assuming that a point zq G dE does not belong to supp/?, hence the 
potential U@ is harmonic in a small disk centered at zq, we arrive to a contradiction by comparing 
the resulting harmonic extension of U 13 with the one given in (16.140 . 

In view of the connectedness of the complement of E and the fact that the support of [3 is 
contained in E the equality U a {z) = U^(z), for z G C \ E, is immediate from Theorem 16.11 (iv). 
Hence supper C E. Furthermore, since the boundary of the domain C \ E in the fine topology 
coincides with its boundary in the Euclidean topology (see e.g. [231 Cor. 1.5.6]), we conclude that 
the equality between the potentials persists in C \ E. 

The last assertion in the corollary can be deduced from Theorem 16.11 (iv)-(v), because this 
guarantees the existence of a cluster point a of the sequence vp n such that U u = U@ on both sides 
of Y\. More precisely, U a = U@ in V \ T\, where V is a neighborhood of Gi not meeting the other 
islands, and therefore a = (3 in such a neighborhood. Similarly we argue for L 2 j_- □ 



8.7. The lemniscate example. 
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Proof of Lemma 17.11 Let (■yw + 1) T / 2 denote the analytic branch in ID = {w : \w\ < 1} that 
equals 1 at w = 0. Then applying Green's formula we have, for j = 0, 1, . . . ,n — 1, 

n fa i \1 lTv 7 dA(w) f j3 n {w) rr\3-r/2 , > 

= y /3 W H( 7 ^ + 1) p^yp = J - (tv + 1) dA(w) 

O D 

7 |7W + 1| T J |7^ + 1| T 

\w\=l \w\=l 

where we have ignored nonzero constants, and in the last equality, we used that (jw + 1) = (j/w+1) 
for \w\ = 1. Consequently, (3 n (w)(j + w) is a monic polynomial of degree n + 1 that vanishes at 
to = —7 and is orthogonal to all polynomials of degree less than n with respect to |du>|/|7u> + l| r . 
The same is true of the right-hand side of (|7.8p and hence the difference of these two polynomials 
(which is of degree < n) must be a multiple of t n (w) that vanishes at —7. Since t n (— 7) 7^ 0, the 
difference of the left and right-hand sides of (|7.8p must be identically zero. □ 

Remark 8.1. It is essential that the cases r = 2,4, ... ,2n be excluded in Lemma 17.11 Indeed for 
r = 2j, where j is a positive integer, it is well-known (cf. [H], §11.2) that t n (w) = w n ~ J (w + 7) J 
for n > j, so that t n (— 7) = in this case. There appears, however, to be no simple formulaQ for 
the polynomials (3 n {w) for such values of r. We shall show in Lemma 18.41 that if r is not an even 
integer, then t n (—^) 7^ for all n sufficiently large. 



Proof of Proposition 17.21 . Here we use the minimality property of the monic Bergman polyno- 
mials Pkm+s( z ) = zS Qk,s( zm )- More precisely, q^^ solves the extremal problem 

I kjS := min{ J \z s q{z m )\ 2 dA : q(t) = t k + ■ ■ ■ G V k }. (8.34) 

Clearly, 



G 



\z s q(z m )\ z dA = m J \z s q(z m )\ 2 dA, 

G Gm 

and the change of variables w = (z m — l)/r m , which maps G m conformally onto the unit disk D in 
the u>-plane, yields 

1 s , m M2 JAI \ r 2m f \q(r m w + l)\ 2 ... . 
\z s q{z m )\ 2 dA{z) = / ' ' , ,w dA (w), 



where 



2 2s . 

r:=2 . (8.35) 

m m 



^For the weight dA/\-fw + 1| 2 , we have 

/3i(w) = w + - + 1- 

7 ln(l-|7| 2 ) 
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Consequently, 

Jim f \ n ( r m„.. I 1 \ 1 2 

I k> „ = — min{ / m 2-,,;' dA(w) : g(t) = t k + ■ ■ • 6 P fc }, (8.36) 
m 7 |r m t(j + l| T 

D 

and, moreover, r~ rnk qi ts (r m w + 1) is the monic (in w) orthogonal polynomial with respect to the 
weight dA(w)/\r m w + l\ T on D. Applying Lemma [7.1l then yields formulas (|7.5p and (|7.6p . provided 
that 7Tfc iS (— r m ) is not zero. In the next lemma we show that this condition is indeed satisfied for k 
sufficiently large. □ 

Lemma 8.4. Let -k^ s (w) be as in Proposition \7.S\ and r be given by 118. 35\) . Then, for each 
s = 0, 1, . . . , m — 2, we have 



,W 2 



r" 



l_/r\ b„ _/l 



7T 



i-^^A-n = Mm/2) ^F(-) +j + 0[ ¥ ) (8.37) 



as k — > oo, where b s is a constant independent ofk. 



Proof. As in [14] . we utilize the results of [12] for Szego polynomials with respect to an analytic 
weight on \w\ = 1. For the weight \w + r m \~ T = l/\r m w + l| r , we have, imitating the notation 
of [12], the following formulas for the exterior and interior Szego functions D e ^ T (w) and Di kT (w), 
respectively, 

D e!T (w) = (^—) , D itT (w) = 0. + r m w)- T ' 2 , (8.38) 



w 

where the branches of the square roots are chosen so that D e ^ T (co) = Di T (0) = 1. The scattering 
function S T (w) is given by 

SJw) = D eT (w)D iT (w) = ( W + rm ) T/2 (l + r m w)- T / 2 for r m < \w\ <r~ m . (8.39) 

V w / 

As shown in [12] (see Equations (16), (25), and (39)), we have for \w\ < rj, where r m < r\ < 1, 

1 r t k S (t) 

D i>T (w)-K k , s {w) = — & ry 1 dt + 0(n 3k ), as k -> oo. (8.40) 

ZTTl J t — W 
1*1=1 

For u; = — r m , we can deform the unit circle in the integral in (18.4Q[) so that the integration takes 
place along each side of the branch cut of D e T (w) joining — r m to to obtain 



J t + r m 

|*|=1 \[-r"\0] [0,-r m ] 



where we utilize the limiting values from below for S T in integrating from — r rn to and the limiting 
values of S T from above in integrating from to — r m . Thus we get (cf. (j8.39|) ) 

7 fc = 2tsin(T7r/2) y | x | T / 2(x + rm) i- r / 2 ^> 
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and on making the change of variable x = — r m (l + cos 9)/2 we find that 



|^ S in(W2)(-l) fc | e- k ^q(0)d9, 



where p(9) := — log(l + cos#) and 

Jim 



1 



-(1 + COS0) 



-T/2, 



'1 + COS I 



\1— T/}T— 1 



sin6'\ r - 1 
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?.42) 



?.43) 



We now apply Laplace's method to deduce the asymptotic behavior of the integral in (|8.42p . 
Since 

oo ^ 

p{9) = - log 2 + Pj el+2 = - log 2 + -0 2 + • • • 

j=0 

and 

oo 

q (0) = Y^q j 9 j+T ~ 1 = (1 - r 2m )~ T / 2 2 1 ~ T 9 T ~ 1 + q 2 9 T+1 + ■■■ , 
j=0 

(note that q\ = 0) we obtain from [17\ Ch. 3, Thm 8.1], that, as k — > oo, 



-kp(tf), 



(0)d0 



F, 2 



t \ (1 - r 2m ) 



2m\-r/2 



+ 



a 2 ' 



fcr/2+1 



+ o 



1 



^r/2+2 



5.44) 



where 02, T is a constant independent of k. From (j8.40p48.44p (taking r] such that < r m < rj) we 
deduce (jOTjh □ 
As an immediate consequence of the preceding lemma we obtain that 



7rfe, s (-r- m ) 



1-7^ + ^ 



2 k 



k 2 



as k — > 00. 



5.45) 



Proof of Proposition 17.31 For s = m — 1 the assertion is obvious from (|7.5p . For |z m — 1| > r 2m 

and s = 0, 1, . . . , m — 2, we appeal to the well-known fact regarding exterior asymptotics of Szego 
polynomials (see e.g. |12j . Proposition 1) that for > r m we have 



lim 

fc^oo W K 



D P Jw) 



w + r m. T /2 



W 



(8.46) 



where the convergence is locally uniform and takes place with a geometric rate. Thus from (|8.45p 
and the representation (|7.6p we deduce (|7.10p for \z m — 1| > r 2m . 

For \z m — 1| < r 2m , we begin with the asymptotic analysis of itf; s (w), for s = 0, l,...,m — 2 and 
\w\ < r m . Assume at first that w ^ [— r m ,0], and consider the integral in the representation (|8.40p . 
For each e > sufficiently small, we can write 



Jk(w) 



1 

27Ti 



t — w 2iri 



1*1=1 



/ 



\t-w\ 



+ + 



[-r m ,0] [0,-r m ] 



t k S T (t) 
t — W 



dt, 



5.47) 
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where integration along both sides of the branch cut from — r m to is as in the proof of Lemma 
18.41 From Cauchy's formula and the representation of S T (t) along each side of the branch cut, we 
deduce that 



J k (w) = w k S T (w) + -sin(r7r/2) 

7T 



m.\r/2 



x k {l + r m x)- T l\x + r m ) 
|x| T / 2 (:r — w) 

which, upon performing the change of variable x = — r m (l + cos#)/2, yields 



■dx, 



i r m(fc+l) 

J k (w) = w k S T (w) + -sin(rvr/2)(-l) fc+1 — j—- 



where p(9) = — log(l + cos 9) and 

m ■■= 

Since 



[i _ r^(i + cos e)]- T / 2 (^-y+ l 9 T+l 

[^(1 + cos #)+«;] (1+ cos fl) T 



{l _ r 2m r r/2 



q(9)=E^ +(T+2yi 

(r rn + w)l T 

(note that q\ = 0), Laplace's method yields 



2 k 



rr(|)(i 



„2m\-r/2 



1 



r m _|_ yj 



+ g 3 T+3 + • • • 



1 £,2+r/2 



£,3+t/2 



3.48) 



as — > oo, where b s {w) is a constant independent of k. Thus, from (|8.48p and (|8.40p . we obtain 
k 1+T / 2 (-l) k+1 sin(r7r/2)rr(T/2) 



Di, T {w)ir k ,s(. w )- 



1 + *£Mf 



r m(fe+i) 27r(l — r 2m ) T / 2 (w + r m ) L 

as A; — ► oo, provided |w| < r m and rf < r m < rj, while for \w\ = r m , w ^ —r m , we obtain 

D i,T{W) L -r = S r (w) + O 



ur 



k l+r/2 



.49) 



3.50) 



as k — > oo, where we take r' m < 77 < 1. 

Combining (|8.45|) with (|8.49|) and (|8.5U|) . we deduce from the representation (|7.6|) that (|7.11|) 
holds for \z m - 1| < r 2m , z m f [1 - r 2m , 1], and that (^T0]) holds for \z m - 1| = r 2m , except for the 
m roots (1 - r 2m ) l l m . In deriving ([7XT]) we used the fact that (z m ) T /V = ^-l^Xs+^/m for 
z G Gj (recall (17. ip ). Finally, by a slight modification of the above analysis, it is easy to see that 
(|8.49p is valid also for w S (— r m , 0] and so (|7.1ip holds for all z satisfying \z m — 1| < r 2m . □ 



Proof of Proposition 17.11 We use the obvious fact that 



km+s 



\Pkm+s{z)\ 2 dA{z). 



3.51) 



G 
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For s = m — 1, we have from (I7.5p . 

r 2 

km+m— 1 



AZ* ! = / Iz™" 1 ^" 1 " = m / I^Cz" - l) k \ 2 dA{z) 



G Gm- 

r 2m(fc+l) r 7rr 2m(fc+l) 



m J m(k + 1) ' 

D 

where, as in the proof of Lemma [8.41 we have made the change of variables w = (z m — l)/r m . Thus 



, m(k + 1) 

Afcm+m-i = y vrr2m(fc+1) • (8-52) 

Now suppose that < s < m — 1. Then, on utilizing the formula (|7.6p we deduce that, for k 
sufficiently large, 

„2m r \ n (^rn„,, _i_ i\|2 



y |z )| dA(z) = — y [rm ^ + 1 | T 



G,7 



2m(fc+l) f \n k+l (w)- ^}[-p 7r k (w)\ 2 



nM^-^—dA( W ), (8.53) 



m J \w + r m | 2 |r m u' + l| r 

D 

where for simplicity of notation we have written ir k = iTk,s- On using the orthogonality property of 
the 7iVs we can simplify the last integral in (|8.53p to obtain 

_ 2 -7r k+1 (-r m )r 2mk+ ™ f \7T k ( W )\ 2 



\m+s 7rfc (_ r m) 2m (A; - 1 + 1) J \r m W + 1^^' 

\w\ = l 

Finally we note that the integral on the right-hand side of (|8.54p equals /J- k 2 , where pt kjS is the 
leading coefficient of the orthonormal polynomial with respect to the weight \dw\/\r m w + l| r on 
the unit circle. As is well-known (see e.g. j!2j . Corollary 2) 

^1,8-^1 = 0(v 2h ) asfc-oo, 
where r m < rj < 1. Combining this fact with (|8.54p and (|8.45p yields (|7.4p . □ 
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